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EXECUTIVE SUMMARY 


ES.l INTRODUCTION 


E3.1.1 Background and Objectives 

A free piston Stirling engine (FPSE) executes a practical 
Stirling cycle without having a crank mechanism. Instead, the 
motions of the two reciprocating elements, the piston and the 
displacer are intimately related to the thermodynamics and 
fluid dynamics of the engine working gas. 

The FPSE has many potential advantages over conventional 
kinematic Stirling engines. They are: 

• Simplicity 

• Freedom from leakage 

• Potential self starting 

• Low noise 

• Long life. 

A hydraulic power system powered by a FFSE may be a very 
attractive alternative for stationary or vehicular applications. 

Sunpower, Inc., of Athens, Ohio had designed and built a 
1 kW FPSE (RE-iOOO) for research testing by NASA-Lewis. As 
shown in Figure ES-1, the RE-1000 had a dashpot in order to dis- 
sipate the output power. 

The objective of the present program is to design and build 
a suitable hydraulic output mechanism for the RE-1000 FPSE to 
demonstrate feasibility, efficiency, and potential practicality. 
This report describes the first phase of the program, design of 
modifications to this particular engine. 


ES.l. 2 Design Specifications 

The present design of the hydraulic output mechanism has 
to satisfy the following requirements: 

• The design should be easily integrable with the 
government-owned RE- 1000 FPSE. 

• Preferably minimal or no modification of the hot section 
of the RE-1000 should be required. 



ORIGINAL PAGE 13 
OF POOR QUALITY 


POWER PISTON 


COOLER 


DISPLACER ROD 


DISPLACER 


REGENERATOR 


HEATER 


Figure ES-1. - Cutaway view of the RE-1000 Free Piston 
Free Displacer Stirling engine built by Sunpower. 


• The RE-1000 FPSE should be coupled to the hydraulic 
fluid through a flexible hermetic interface. 

• The following pressure and temperatures are to be used 
for design point. 

1. neater tube temperature, 650°C 

2. Coolant temperature, 30°C 

3. Working gas charge pressure, 70 bar. 

• The design approach should be scaleable over the power 
range of 65 to 400 kW. 

• The flexible interface should have a fatigue life in 
excess of 500 hr of operation for test purposes; yet 
have the potential to design for much longer lifetime. 

ES.1.3 Program Description 

The program consisted of three tasks. Table ES-1 shows 
each of the tasks divided into specific subtasks. 


TABLE ES-1. - DESCRIPTION OF THE PROGRAM EY TASKS 


Task 

I 

Preliminary Analysis and Concept Selection 


1.1 

Preliminary Systems Analysis 


1.2 

Analysis and Design of Diaphragm/Bellows 


1.3 

Development of Three Candidate Design Concepts 


1.4 

Evaluation and Comparison 


1.5 

Design Review Meeting 

Task 

II 

Detailed Design 


II. 1 

Detailed Design Analysis 


II. 2 

Design Optimization 


II. 3 

Sensitivity Analysis 


II. 4 

Detail Design Review Meeting 


II. 5 

Final Design and Detail Drawings 

Task 

III 

Reporting Requirements 


III. I 

Monthly Reports 


III. 2 

Design Review Meetings I and II 


III. 3 

Final Report 



ES.2 PRELIMINARY DESIGN 


ES.2.1 Description of System Concepts 

Three candidate system concepts were prepared. Each are 
described below. 

Candid a te design No. 1. - As shown in Figure ES-2, this 
design satisfies the requirements and objectives of the program. 
This design requires a single diaphragm, two accumulators, and 
two check valves. One disadvantage of this concept is that 
the maximum hydraulic pressure is limited to approximately the 
maximum engine pressure. 

Candidate design No. 2. - This is the simplest design with 
diaphragm/bellows interface that can provide the high pressure 
hydraulic output directly as shown in Figure ES-3 without using 
an inefficient hydraulic pump-motor combination. 

Candidate design No. 3. - The design shown in Figure ES-4 
is the best near-term candidate. It does not involve any major 
changes in the RE-1000 FPSE and even retains the solid piston. 
The gas isolation bladder seems much easier to design than dia- 
phragms or bellows for design No. 1 or No. 2 since there will 
be pract' rally no pressure differential across the bladder. 


ES.2. 2 Dynamic Stability and the Pump Design 

An important consideration in the design of a hydraulic 
output mechanism for a FPSE is the dynamic stability of the com- 
bined system. Figure ES-5 shows a power load matching between 
an ordinary hydraulic pump and the RE-1000 F°S£ assuming a con- 
stant frequency. Because the slope of the FPSE power curve at 
the matching point is higher than that of the hydraulic pump 
load, any deviation of the piston stroke from the matching 
point will be amplified by the surplus in power if the stroxe 
should increase, or by the shortage in power if the stroke should 
decrease . 

Consequently, the system consisting of the RE-1000 F> SE 
and an ordinary hydraulic pump is inherently unstable. 

In order for a hydraulic output to be inherently stable, the 
power load curve near the match point should look like the one 
shown in Figure ES-6 where the slope of load curve is higher than 
the slope of the engine power curve. In this case, any deviation 
of stroke from the matching point will be corrected by the shortage 
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Figure ES-3. - Candidate design No. 2 with a diaphragm 
interface transformer piston, and a null-center 
band pump output. 
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Figure ES-4. - Candidate design No. 3 with the engine 
piston, an isolation bladder, and a null-center 
band pump output. 
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Figure FS-5. - Power versus stroke curves of a FPSE and an 
ordinary hydraulic pump showing an unstable 
operating point. 



Figure ES-6. - Power versus stroke curves of a FPSE and an 
acceptable hydraulic pump showing a stable 
operating point. 
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in power if the stroke should increase, or by the surplus in 
power if the stroke should decrease. 

The three hydraulic output methods shewn in Figures ES-7, 
ES-8, and ES-9 satisfy the above criteria and are inherently 
stable. The design shown in Figure ES-7 is for the candidate 
design No. 1, whereas the designs shown in Figures ES-8 and ES-9 
are applicable to both the candidate designs No. 2 and No. 3. 

The volume displacement of the stabilizer diaphragm shown 
in Figure ES-7 is less waan that of the main engine diaphragm. 

If the main diaphragm displacement is less than that of the 
stabilizer diagram, there will be no pumping. In other words, 
this pump has a dead band in the intermediate stroke region. 

Figure ES-8 shows the schematic of a "double" acting 
hydraulic pump which also has a mid-stroke hydraulic dead band. 
The dead band is achieved by properly positioning port holes. 

Figure ES-9 shows the schematic of a "single" acting pump. 
It uses essentially the same concept as the "double" acting pump 
except that it has a mean pressure reservoir. 


ES.2.3 Comparison of System Concepts 
and Final Selection 

Of the three candidate design concepts described in the pre- 
vious section, the candidate design No. 1 was eliminated. The 
reasons for its elimination are: 

• Limited pressure output. 

• Low efficiency if a hydraulic motor pump is used for 
pressure boosting. 

The remaining concepts to be examined were candidate design 
No. 2 having a flexible interface with the working gas (diaphragm 
or bellows) , and the candidate design No. 3 with the piston 
retained. Table ES-2 presents a qualitative comparison of the 
concepts by listing their relative advantages and disadvantages. 

The bellows version was discarded because of the following 
reasons : 

• High projected cost in manufacturing and quality 
control . 
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Figure ES-7. - Liquid piston driven hydraulic pump with a 
null- center band. 
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Figure ES-9. - Single acting hydraulic pump with a null' 
center band. 


TA3LE ES-2. - COMPARISONS OF CANDIDATE DESIGNS OF THE HYDRAULIC 
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• Difficulty of providing overpressure protection. 

• Fluid squeezing loss problems. 

Of the two remaining candidate designs, the diaphragm ver- 
sion was selected to be developed further. The rationale was 
that the diaphragm version would provide more opportunities 
for technical advancement than the piston version. 


ES.2.4 Preliminary Dynamic and Thermodynamic Analysis 

The objectives of the preliminary analyses were: 

• To check the compatibility of the Fester-Miller designed 
hydraulic pumps with the RE-1000 FPSE and to determine 
if the overall system is dynamically stable. 

• To obtain the performance trend for the RE- 1000 FPSE 
integrated with the pump. 

• To help the design effort by indicating necessary 
modifications as the design evolved. 


ES.2.4.1 Dynamic analys is. - A FPSE dynamic analysis, which had 
been developed by Sunpower, was used to investigate the dynamic 
compatibility of the RE-1000 FPSE to particular hydraulic loads. 
The analysis utilizes simplified thermodynamics. 

From the dynamic analysis, the following were learned: 

• The system consisting ot the RE-1000 FPSE and Foster- 
Miller designed hydraulic pump is dynamically stable. 

• At nigh power output, the displacer collides with tne 
piston. This can be prevented by shortening the piston 
or moving its position. 

• The nonsinusoidal motion of the piston can be corrected 
by increased stored energy in the system. 

• The analysis was for the piston version of the RE-1000 
FPSE, but the trends should also apply to the flexible 
interface version with equivalent inertial energy. 


ES.2.4. 2 Thermodynamic ana l ysis . - A number of computer simula- 
tions were run of the RE-1000 FPSE engine to determine the effects 
of increased heat transfer in the compression space resulting from 
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the substitution of a -diaphragm for the current piston. The 
simulations employed Sunpower’s third order code with constrained 
dynamics of piston and displacer. 

The conclusions of the preliminary thermodynamic analysis 
were as follows: 

• All components of the Stirling cycle must be considered 
together before any conclusions can be drawn regarding 
the merits of an increased heat transfer to approach 
isothermalization . 

• Merely increasing the heat transfer in the compression 
space of a Stirling engine is not sufficient to guaran- 
tee an improvement in performance. The machine must 

be reoptimized in all its aspects before near isother- 
malization can be expected to yield a benefit. 

• The heat transfer in the diaphragm area should be 
carefully investigated, because design and performance 
depend strongly on prediction of this phenomenon. 

• Net heat transfer from, or cooling of the diaphragm 
compression space appears necessary to achieve good 
performance . 
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ES.3 DETAIL DESIGN 


ES.3.1 Description of the Recommended Design 

The final design incorporates the following: 

• Modular sections and simple, accessible parts for 
research testing. 

• Interchangeability with current RE-1000 FPSE output 
section without disturbing engine section. 

• Either vertical axis hydraulic assembly, or horizontal 
assembly to demonstrate dynamic balancing better. 

• Displacement sensors on all engine moving parts and 
hydraulic pump. 

• Conservative design for low stress and low cost. 

The design is described referring to Figure ES-10, which is a 
simplified schematic of the actual design. Basically the system 
operates as follows: 

• Engine gas displaces oil trapped between the engine 
diaphragm and the bottom end of the transmission piston. 

• The transmission piston displaces oil trapped on the 
other (top) end, thereby displacing the bounce gas 
diaphragm. 

• No significant net power is transferred to this oil or 
the diaphragms, and the effective mass of this piston 
and oil is insufficient to provide the desired cyclic 
motion in combination with engine and bounce gas springs. 

• Additional effective piston mass is provided in the 
counter balance assembly. 

• The annular piston masses are forced by the surrounding 
trapped oil to move in opposite directions synchronously 
so that the entire engine system is effectively balanced 
dynamically. 

• If only one mass is clamped to the piston rod and oil 
is free to flow around it in the same section, then 
unbalanced dynamic operation results. 
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Figure ES-10. - Simplified schematic of the recommended 
hydraulic output system design. 
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The piston rod transmits net power to the hydraulic 
oil pump, which is isolated from the diaphragm oil. 


It is obvious that the basic layout affords easy access for 
modification. 


ES.3.2 Detailed Analysis 

The detailed analyses include the metal diaphragm analysis, 
the thermodynamic and dynamic analyses of the engine, and 
estimates of additional power losses due to hydraulic output 
conversion of the RE-1000 FPSE . 


ES .3.2.1 Metal diaphragm analysis . - Two metal diaphragms are 
used in the hydraulic output version of RE-1000 FPSE. They are 
the engine diaphragm and the bounce gas diaphragm. The engine 
diaphragm provides a hermetically sealed interface between the 
engine gas and the hydraulic fluid. The bounce gas diaphragm 
separates the hydraulic fluid from the bounce space gas. 

The diaphragm analysis provided the information on the 
following items of interest: 

• btress level and distribution. 

• Axial displacement distribution. 

• Meridional displacements of edges. 

• Stiffness of diaphragm near TDC and BDC . 

For the analysis of the diaphragms, two existing computer 
programs, NONLIN ind ANSYS, were utilized. The program "NONLIN" 
is a bellows/diaphragm ^analysis developed by Dr. J.F. Lestingi, 
et al. (ref. 1) at Battelle Memorial Institute, Columbus, Ohio. 

The "ANSYS" program is a proprietary general purpose finite 
element program developed by Swanson Analysis, Inc. 

Table ES-3 shows the specifications and boundary conditions 
of the two diaphragms used in the analysis. 

Figure ES-11 shows the stress versus cycles to failure data 
for the typically recommended AM 350 SS sheet of 1.016 mm in 
thickness. The maximum stresses of the diaphragms are 100 'v 
110 MPa which gives the safety factor of at least five with 
respect to the endurance limit of 590 MPa for AM 350 SS. 
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Engine Diaphragm 


Bounce 

gas diaphragm 


Specifications 



Outer diameter 

(cm) 

mm 

20.32 

Inner diametec 

(cm) 

6.35 

2.45 

Thickness 

(cm) 

0.0508 

0.0508 

Youngs Modulus, E 

(MPa) 

2.069 x 10 5 

2.069 x 10 5 

Poisson’s Ratio, v 

0.27 

0.27 

Boundary conditions 

Outer edge 

Simply supported 
to the casing 

Simply supported 
to the casing 

Inner edge 

Simply supported 
to the center 
plug which moves 
with the 
diaphragm 

Axially clamped 
to and radially 
guided by the 
center plug which 
moves with the 
diaphragm 

Volume amplitude 
at design point 
(cm 3 ) 

40 

40 





























800 


ORIGINAL PAGE IS 
OF POOR QUALITY 



♦ 






-J 






ft 

<o 

o. 

3: 



/I 

O 

rH 

O 


o 


< x 

/ f 

e— « 


LO 


■ 

/ 



ro 


= *— O 

/ * 


• 

rH 


Q.ZM 



CD 




' 3 


Z 

II 


«OK / 

LU <£ 


»— « 



Qa< / 

U Ct Z 


u 

ZD 



ZlOh 



h- 


LU 2: W 

<t Q. —• Z 


0 

Li. 


z ua 

_J «£ OO •— > 


ae 


CD 

< — * » — « >_ 

<C •— 1 LU O 



•» 

z 

CD 00 

CQ Q Q Q. 


Li. 

H 

►H 

Z LlI 0 


Q 

O 

O LU 

Q — . 

UJQH 


r -H 


U") LU 

Z ^ 




Z 

CO z 

LU 




O 

l/) 

CD LU 




»— « 

H 

O 





O E 

o z 




O 

co E 

LU LU 




LU 


IS) cl 




a: 

O vo 

CL LU 




►H 

*<) f-H 

LU Ll. 




O 

CO O 

> LU 





3; • 

LU CL 





< •—« 

CL 




2 






z 





O 

►— * 
0 

ZD 


l/> 


CJ 

>- 

o 


<o 

o 


o 

z 

o 


CD 

z 

LU 

ae: 

GO 


in 

c 


h- 
*— « 



4J 

0) 

a> 

x: 

m 

U) 

V) 

o 

IT) 

m 


M 

0 


CJ 

u 


•H 

ifl 

y-i 


o 

(0 

0) 

o 

>1 

0 

U) 

3 

(0 

d) 

> 

a 

w 

( 1 ) 

M 

■P 

W 

1 


I 

CO 

w 

<u 

u 

3 

CP 

*H 

Cl4 


PdW SS381S 


20 


This means that as long as the diaphragms are fabricated without 
defects, they will last indefinitely under normal operating con- 
ditions. This probably is also true for many other more con- 
venient ..heet materials. 


ES . 3.2.2 Thermodynamic performance predictions . - For the detailed 
analysis of the hydraulic output Free Displacer Stirling engine, 
the Sunpower Third Order Stirling engine analysis was utilized. 

This has two modes of operation: the constrained mode uses pre- 

determined dynamics of reciprocating components. The uncon- 
strained mode simultaneously solves the dynamics and thermo- 
dynamics. The highly impulsive or step-function-l'ike nature of 
the load necessitated the use of the unconstrained mode to closely 
simulate the actual processes of the engine. The less expensive, 
constrained mode analysis was used to obtain the performance trend 
and to perform a sensitivity analysis of important operating 
parameters . 

Within the scope of the design proonm, there are two para- 
meters significantly affecting the power output and the effici- 
ency of the engine: the displacer phase angle and the cylinder 

or diaphragm-space heat transfer. 

The displacer phase angle is closely linked to the overall 
engine dynamics. Varying the displacer phase angle will change 
the frequency slightly and the amplitudes of the piston and the 
displacer, the pressure wave phase angle, and consequently the 
engine output performance. 

With the increased heat transfer area and the extremely 
narrow spaces in the diaphragm region of the compression space, 
estimation v/t the magnitude of the heat transfer has significant 
consequences in the predicted performance of the engine. 

Hence, the analyses were performed iteratively with par- 
ticular attention to the combined effect of displacer phase angle 
and the compression space heat transfer. Note that the diaphragm 
heat transfer faalysis is original, based on limited related 
data; actual levels remain to be demonstrated as a significant 
part of the entire FPSE research program. 

The results of the third order analysis included: 

• The base-point, or reference simulations of the standard 
RE-1000 FPSE and the diaphragm piston engine. 

• The sensitivity analysis of the diaphragm piston 
Stirling engine with respect to the cylinder heat 
transfer, h, the phase angle of the displacer, 4> d , 
and the pump load force, F i oac j* 
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• The predicted performance of the diaphragm piston 
Stirling engine. 

Figure ES-12 illustrates the parametric effect of heat 
transfer in the compression space. Included in Figure ES-12 are 
P-V work, indicated efficiency and availability loss for the 
compression space as a function of the rate of heat transfer. 

The following observations were made: 

• The P-V power near the adiabatic extreme (small hA) 
is higher than that of the isothermal extreme (large 
hA) This is because the mass flow in and out of the 
compression space for the near adiabatic case is 
higher than that of the near isothermal. 

• The efficiency for the near isothermal case is much 
higher than that of the near adiabatic case. 

• In the intermediate region between the adiabatic *nd 
the isothermal case, the efficiency and the powei 
reduce to their respective minimum values. Unfor- 
tunately, the diaphragm engine appears to fall in this 
broad region, and there seems no obvious, simple way to 
avoid it. 

Figure ES-13 shows the efficiency and power of the diaphragm 
engine with respect to the displacer phase angle based on con- 
strained mode third order analysis. The maximum efficiency 
occurs near 60 -deg and the maximum power occurs near 80 deg. 
Figure ES-14 shows a simple modification of the displacer rod 
which would increase the gas spring volume to achieve 8( n 
phase angle. 

In Table ES-4, the predicted engine performances are _ 
sented. They include the results from separate analyses of 
losses due to viscous friction, leakage, and pressure drop in 
the hydraulic section of the engine. The improvement of per- 
formance due to the simple displacer rod modification is notable. 
It is predicted that the hydraulic output Stirling engine, as 
designed, is capable of delivering 1.87 kW brake power at 
29.7 percent efficiency. 

Considering that this particular engine design is not yet 
optimized, the predicted level of performance looks promising. 
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Figure ES-12. - Predicted performance of the Diaphragm Piston, Free Displacer 
Stirling engine as a function of heat transfer in the compression space. 
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Figure ES-14. - Proposed displacer rod design 
increase the gas spring volume. 
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TABLE ES-4. - PERFORMANCE PREDICTIONS OF DIAPHRAGM PISTON 
FREE DISPLACER STIRLING ENGINE WITH NULL 
CENTER BAND HYDRAULIC PUMP OUTPUT 



Displacer phase angle 
(degrees) 


Indicated power 
(watts) 


Indicated efficiency 
(percent) 


Power loss due to 
hydraulic system 
(watts) 

With balance 
mechanism 

Without balance 
mechanism 


Brake power (watts) 

With balance 
mechanism 

Without balance 
mechanism 


Brake Efficiency 
(percent) 

With dynamic 
oalance 

Without dynamic 
balance 


With the high phase 
displacer rod 

With the new 
recommended 
displacer rod 

98 

77 

1210 

2010 

24.2 

32.0 

212 

212 

(142) 

(142) 

998 

1798 

(1068) 

(1868) 

20 

28.6 

21.4 

29.7 






1 . INTRODUCT ON 


1.1 Background and Objectives 

A free piston Stirling engine (FPSE) is a machine which ex- 
ecutes a practical Stirling cycle without a crank mechanism. 

It was invented in the late 1950 's by William Beale. The motion 
of the two reciprocating elements, the piston and the displacer 
is intimately related to the thermodynamics and fluid dynamics 
of the engine working gas. 

The FPSE has many potential advantages over conventional 
kinematic Stirling engines. They are: 

• Simplicity 

• Freedom from leakage 

• Potential self starting 

• Low noise 

• Long life. 

A FPSE can provide electrical output by a linear alternator 
arrangement. It can also provide a hydraulic output by a hydrau- 
lic pump arrangement. A hydraulic power system coupled with a 
FPSE is a very attractive alternative for stationary or vehicular 
applications . 

Sunpower, Inc., of Athens, Ohio, had designed and built a 
1 kW FPSE having a load simulated by internal gas pumping the 
RE-1000 for research purposes at NASA- Lewis. Figure 1 shows the 
cutaway view of the RE-1000. The objective of the present pro- 
gram is to design a suitable hydraulic output mechanism for the 
RE-1000 to demonstrate analytically the feasibility of the hydrau- 
lic output FPSE. 


1.2 Design Specifications 

The present design of the hydraulic output mechanism has to 
satisfy the following requirements: 

• The design should be easily integrable with the 
government-owned RE-1000 FPSE. 

• Preferably minimal or no modification of the hot sec- 
tion of the RE-1000 should be required. 
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Figure 1. - Cutaway view of RE-1000 Free Piston 
Displacer Stirling engine built by Sunpower. 













• The RE-1000 FPSE should be coupled to the hydraulic 
fluid through a flexible hermetic interface. 

• The following pressure and temperatures are to be used 
for design point. 

1. Heater tube temperature, 650°C 

2. Coolant temperature, 30°C 

3. Working gas charge pressure, 70 bar. 

• The design should be scalable over the power range of 
65 kW to 400 kW. 

• The flexible interface should have a fatigue life in 
excess of 500 hr of operation for test purposes; yet 
have the potential to design for much longer lifetimes. 


1.3 Program Description 

The program consisted of three tasks. Table 1 shows each 
of the tasks divided into specific subtasks. 

It is expected that a subsequent program will involve fabri- 
cation and initial hydraulic testing of components to be instal- 
led and tested on the RE- 1000 at NASA-Lewis. 
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TABLE 1. - DESCRIPTION OF THE PROGRAM BY TASKS 


Task I 

Preliminary Analysis and Concept Selection 

1.1 

Preliminary Systems Analysis 

1.2 

Analysis and Design of Diaphragm/Bellows 

1.3 

Development of Three Candidate Design Concepts 

1.4 

Evaluation and Comparison 

1.5 

Design Review Meeting 

Task II 

Detailed Design 

II. 1 

Detailed Design Analysis 

II. 2 

Design Optimization 

II. 3 

Sensitivity Analysis 

II. 4 
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2. PRELIMINARY DESIGN 


In this section, the candidate design concepts will be 
briefly described and be followed by the discussions of the system 
stability, pump design, trade-off studies, and the preliminary 
analyses 


2.1 Description of System Concepts 

Three system concepts are prepared. Each one will be de- 
scribed below. 

Candidate design No. 1. - An shown in Figure 2, this design 
is the simplest one that satisfies the requirements and objectives 
of the program. This design requires a single diaphragm, two 
accumulators, and two check valves. One disadvantage of this 
concept is that the maximum hydraulic pressure is limited to 
approximately the maximum engine pressure. To overcome this dis- 
advantage, an additional hydraulic pump-motor combination can be 
included; however, this will somewhat lower the overall efficiency 
of the system. 

There seems tc be no apparent obstacle to scale this design 
up to, say, 50 kW per cylinder. In order to provide a higher 
power output (up to 400 kW) , multi-cylinder arrangements are 
necessary. In such cases, only a single set of high and low pres- 
sure accumulators will be needed, but the method of phasing the 
cylinders has to be devised. 

Candidate design No. 2. - This is the simplest design with 
diaphragm/bellows interface that can provide the high pressure 
hydraulic output directly without the inefficient hydraulic 
pump-motor combination as shown in Figure 3. It is also scalable 
up to a reasonable size (50 kW per cylinder) . In ganging up 
several cylinders together, a two cylinder unit has a definite 
advantage over other multi-cylinder units in that the two-cylinder 
arrangement eliminates the need for bulky bounce gas accumulators 
by using two cylinders hydraulically coupled out of phase by 
180 deg. Other multi-cylinder units will require the means of 
maintaining proper phase angle between the cylinders. 

Candidate design No. 3. - The design shown in Figure 4 is 
the best near-term candidate. It does not involve any major 
changes in the RE-1000 FPSE and even retains the solid piston. 

The gas isolation bladder seems much easier to design than 
diaphragms or bellows for design No. 1 or No. 2 since there will 
be practically no pressure differential across the bladder. Also, 
the high pressure, low volume transmission diaphragms are very 
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BOUNCE SPACE 
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(DISPLACEMENT '100 cc) 
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ACCUMULATORS (2) 
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PIPE LENGTH 


CHECK VALVES (2) 


LOAD 

HYDRAULICS 


Figure 2. - Candidate design No. 1 with a diaphragm interface 
and a liquid piston driven null-center-band pump output. 
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COMMERCIAL 
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Figure 4. - Candidate design No. 3 with the engine piston 
isolation bladder, and a null-center-band pump output. 
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small in size and will cost much less than those for design No. 1 
or No. 2. Compared to the other two candidate designs, this 
design will pose the fewest questions in scaling up among the 
three designs. 

Agair, the two cylinder arrangement will have definite 
advantages as in design No. 2; it will not require a large gas 
spring volume and will have an inherent phasing characteristic 
of the two cylinders if the gas spring volume is common to both. 


2.2 Dynamic Stability and the Pump Design 

An important consideration in the design of a hydraulic 
output mechanism for a FPSE is the dynamic stability of the com- 
bined system. Figure 5 shows a power load matching between an 
ordinary hydraulic pump and the RE-1000 FPSE assuming a constant 
frequency. Because the slope of the FPSE power curve at the 
matching point is higher than that of the hydraulic pump load, 
a small increment in stroke from the matching point will result 
in surplus power in the system which will cause an increase in 
stroke. The process repeats itself so that the stroke increases 
continuously until impact. Obviously, should the stroke decrease 
just a small amount from the matching point, the load absorbs 
more power than the engine puts out. This tends to load up the 
engine and shorten the stroke. The shortened stroke creates 
further power shortage. The process repeats itself so that the 
stroke decreases continuously until the engine stops. 

Consequently, the system consisting of the RE-1000 FPSE 
and an ordinary hydraulic pump is inherently unstable. To make 
it stable requires a power/load feedback control system, which 
seems disadvantageous due to the complexity, cost, and poor re- 
sponse or inefficiency of typical Stirling engine control methods. 

In order for a hydraulic output design to be inherently 
stable, the power load curve near the match point should look 
like the one shown in Figure 6 where the slope of load curve is 
higher than the slope of the engine power curve. In this case, 
any deviation of stroke from the matching point will be corrected 
by the shortage in power output if the stroke should increase 
or by the surplus in power if the stroke should decrease. 

The three hydraulic output methods shown in Figures 7, 8, 
and 9 satisfy the above stability criteria and are inherently 
stable. Figure 10 shows the dynamic characteristics for the de- 
sign of Figure 7. The design shown in Figure 7 is for the candi- 
date design No. 1, whereas the designs shown in Figures 8 and 9 
are applicable to both the candidate designs No. 2 and No. 3. 
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Figure 5. - Power versus stroke curves of a FPSE and an 
ordinary hydraulic pump showing an unstable 
operating point. 



STROKE 


Figure 6. - Power versus stroke curves of a FPSE and an 
acceptable hydraulic pump showing a stable 
operating point. 
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Figure 7. - Liquid piston driven hydraulic pump with 
null center band. 






The volume displacement of the stabiliser diaphragm snown 
in Figure 7 is less than that of the main engine diaphragm. When 
the volume displacement of the main engine diaphragm surpasses 
that of the stabilizer diaphragm, the check valves open, and the 
pump starts pumping. If the main diaphragm displacement is less 
than that of the stabilizer diaphragm, there will be no pumping. 

In other words, this pump has a dead band in the intermediate 
stroke region. 

Figure 10 illustrates how this pump works. At t ■ tj, both 
diaphragms are at top dead center (TDC) . As the main diaphragm 
moves down, the stabilizer diaphragm moves down also, until the 
latter reaches bottom dead center (DDC) at t ■ to. The main dia- 
phragm starts pumping to the high pressure reservoir at t ■ t£ until 
it reaches its BDC at t - t 3 . Retween t - t 3 and t ■ t 4 , both dia- 
phragms move toward TDC until at t ■ t 4 the stabilizer diaphragm 
reaches its TDC first. The main diaphragm starts suction from 
the low pressure reservoir at t - 14 until it reaches its TDC at 
t t 5 . This constitutes a complete pumping cycle. 

The loading function for the liquid piston shown in Figure 10 
contains not only the power absorbing components but also some 
spring components. This will tend to increase the frequency of 
the FPSE. However, the system will have a stable matching point 
like the one shown in Figure 6 . 

Figure 8 shows the schematic of a double acting hydraulic 
pump which also has an intermediate stroke dead band. The dead 
band is achieved by properly positioning port holes so that 
within a certain stroke amplitude, the device will not pump. 

The leakage passages indicated by broken lines are to prevent 
cavitation during the sudden suction in the opposite chamber 
when the piston begins to cover the port hole. 

Figure l > shows the schematic of a single actinq pump. 

It uses essentially the same concept as the double acting pump 
except that it has a mean pressure reservoir. 

Figure 11 shows the load function for the double or single 
actinq pump in relation to the pump piston displacement. For 
t he single acting pump, the load is proportional to (Fjj-Pj )/ 2 . 

For the double acting pump, the load is proportional to P H -P L . 

This means the double acting pump will be much smaller in size 
than the single acting pump for the same capacity. 

In Figure 11, at t 0, the piston is at the midstroke 
position. At, t * tj, the piston covers the port hole and begins 
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Figure 10. - Voluzoe displacement and loading versus time for the liquid piston 
driven hydraulic pump with a null center band. 





to pump to the high pressure (P H ) reservoir until the piston 
reaches TDC at t * tj. As the piston moves away from TDC , it 
draws in from the low pressure (P^) reservoir until at t = tj 
when the port hole is uncovered and the suction process stops. 
Similar processes occur as the piston travels in the opposite 
direction. The loading function shown in Figure 11 does not con- 
tain any spring components as were present in the liquid piston 
case . 


Figure 12 shows the power matching between a FPSE and the 
hydraulic output methods presented in Figures 7, 9, and 10. The 
dead band width is dependent on the location of the port holes 
or the relative magnitude of the stabilizer diaphragm displace- 
ment. The slope of the idead load curve is directly proportional 
to the effective bore area of the pump and the accumulator 
pressure. 

In Figure 13 are presented the various power loss mechanisms 
and their relative magnitudes for the proposed hydraulic pumps. 
The leakage control method to prevent cavitation in the double 
acting pump should not be a significant source of inefficiency. 

In order to design an efficient hydraulic output system, all of 
these losses have to be carefully controlled. 


2.3 Comparison of System Concepts and Final Selection 

Of the three candidate design concepts described in the 
previous section, the candidate design No. 1 was discarded. The 
reasons for its elimination are as follows: 

• Limited pressure output 

• Low efficiency if a hydraulic motor is used. 

The remaining concepts to be examined were candidate design 
No. 2 having a flexible interface with the working gas (diaphragm 
or bellows) , and the candidate design No. 3 with the original 
piston retained. Table 2 presents a qualitative comparison of 
the concepts by list^ g their relative advantages and disadvantages. 
The items of comparison are approximately in descending order of 
importance . 


2.3.1 B ellows version. - A bellows in the present application 
has to satisfy a very difficult set of requirements, which 
include the following: 
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Figure 1 
output 



- Power load matching between a FPSE and the 
hydraulic output system design. 



. - Power losses in the proposed hydraulic 
system design - noc to scale. 



TABLE 2. - COMPARISONS OF CANDIDATE DESIGNS OF THE HYDRAULIC 
OUTPUT SYSTEM FOR RE-1000 FPSE 




























































• High frequency. 

• Low dead volume or clearance volume. 

• Low gas/oil flow losses. 

• Overpressure protection in nested condition as well 
as in fully extended position. 

A leading manufacturer of metal bellows, who also markets bellows 
type compressors, was consulted about the possible use of metal 
bellows in the hydraulic output applications. None of the existing 
bellows designs could meet all of the above requirements. Also, 
it was learned that non-standard bellows might require a new die 
for fabrication, and the weldment geometry and stresses do not 
scale simply. (Formed bellows are impractical for this purpose 
due to large size and/or clearance volume) . 

The bellows version was discarded because of the following 
reasons : 

• High projected cost in manufacturing and quality 

control . 

• Poor scalability. 

• Difficulty of providing overpressure protection. 

• Fluid squeezing loss problem. 

Therefore the final choice has to be either the candidate 
design No. 2 with a diaphragm (diaphragm version) or the candidate 
design No. 3 (piston version) . These two will be discussed 
further . 


2.3.2 Diaphragm version. - The preliminary design of the diaphragm 
version has the following features: 

• Annular diaphragm with simply supported edges; that is 
clamped between O-ring seals. 

• Overpressure protection mechanism provided for both 
directions on both diaphragms. 

• Minimal flow loss due to diaphragm and support 
structures . 

• Active cooling of both gas side and oil side provided 
for engine diaphragm. 
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• Displacement transducers are on two diaphragms and 
engine displacer. 

• Dynamic mass balancing having moderate hydraulic 
losses . 

• Heavy, but low cost experimental version of diaphragm 
engine. 

The diaphragm version will demonstrate four areas of interest 
hydraulic output mechanism, diaphragm performance, mass balance, 
and hermetic gas/oil separation and gas containment. The risks 
involved with the diaphragm version are the relatively high cost 
to convert the RE-1000, potentially poor engine performance, and 
the large test program required to demonstrate system performance. 
The potential for poor performance exists primarily due to the 
generally inadequate state of the art with respect to cyclic 
heat transfer in the diaphragm gas compression space, as discussed 
in subsection 3.2. 


2.3,3 Piston version. - The piston version has the following 
features : 

• Only small amount of modification of RE-1000 required. 
Use existing piston, pressure vessel. 

Gas/oil separation achieved by a simple method of 
using- a pre-stretched rubber tube or drum head clamped 
to piston rod. 

» Rubber failures will not cause any immediate oil con- 
tamination of bounce space. 

• All hydraulic components except piston and rod are 
located outside the hermetically sealed pressure vessel. 
Adaptation of check valves, hoses, etc., do not involve 
depressurization of the engine 

The piston version would demonstrate two areas of interest: 

• Hermetic gas/oil separation and oil containment. 

• Hydraulic output mechanism. 

The risks involved in selecting the piston version are minor 
compared to the risks of the diaphragm version. The risks associ- 
ated with the piston version are related to the gas/oil separation 
bladder design and the small hydraulic diaphragm development. The 
advantages of the piston version are: 
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• Low development and production cost for conversion to 
hydraulic output. 

• Minimal test program required. 

• No change of engine performance (except what is generic 
to hydraulic output) . 

It was obvious that selecting the piston version over the 
diaphragm version would mean the least developmental effort and 
risks. However, both versions could be developed within reason- 
able cost and time. 

It was decided to select the diaphragm version to be developed 
further. The rationale was that the diaphragm version would pro- 
vide more opportunities for technical advancement. For example, 
the diaphragm version development program would provide the 
following : 

• Modularized design for independent development of 
output components. 

• Data on actively cooled, diaphragmed cold-end FPSE. 

• High cyclic rate data on diaphragms. 


2.4 Preliminary Dynamic and Thermodynamic Analysis 

The objectives of the preliminary analyses were: 

• To check the compatibility of the FMA designed hydraulic 
pumps with the RE-1000 and determine if the overall 
system is dynamically stable. 

• To obtain the performance trend for the RE-1000 
integrated with the pump. 

• To help the design effort by indicating necessary 
modifications as the design evolved. 

To achieve the objectives, Sunpower performed both dynamic 
and thermodynamic analyses. The results of the analyses will 
be presented here. 


2.4.1 Dynamic analysis. - This subsection summarizes the results 
of dynamic analysis performed by Sunpower to investigate the 
dynamic compatibility of the RE-1000 to particular hydraulic 
loads . 
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The dynamic studies utilized a Sunpower developed pressure 
factor analysis. This analysis, while not intended to give ex- 
act numerical results, is a useful tool for predicting approxi- 
mate engine dynamics and performance trends. 

The general trends of the pressure factor analysis had been 
verified, in some cases, using the Sunpower third order analysis. 

The RE-1000 now has two sets of displacer and rods. For the 
present analysis, the set designed for maximum power was used in 
determining the hydraulic load-engine compatibility. The RE-1000 
with this set had run with a displacer stroke of 2.90 cm, a phase 
angle of 80 deg, and a frequency of 34 . 7 Hz. The hydraulic pump 
is designed to pump only if the stroke becomes larger than the 
dead band width of 2.748 cm. 

The load function to simulate the pump load is similar to 
the one shown in Figure 11. Table 3 presents the results of 
the analysis. Cases No. 1 through No. 4 show the effect of the 
increasing load on the overall performance. 

Figure 14 illustrates the effect of load force on the per- 
formance of the hydraulic output, RE-1000. The square wave load, 
with a center null band, has a tendency to slow down the motion 
of the piston near the extremes of its travel. This explains 
why the frequency decreased as the load force amplitude increased 
as shown in Figure 14. 

It was aiso learned that the square wave load affects the 
dynamics of the piston and displacer significantly. The strokes 
of the piston and displacer change and their motions are dis- 
torted from sinusoidal. For case No. 1, the displacer collides 
with piston as indicated in Figure 15. 

Because this displacer was designed for near resonant 
operation, decreasing frequency with load reduces the displacer 
and piston amplitudes. This was true for cases No. 2, 3, and 4 
when there were no collisions predicted. 

The displacer-piston collision during high power output can 
be prevented by shortening the length of the piston. The addi- 
tional compression space volume should result in only a small 
reduction in engine power and frequency. 

In an attempt to bring the piston displacer motion closer 
to a sinusoidal condition, the energy stored in the piston was 
increased by increasing the piston mass. Starting with case No. 5 
through No. 8 in Table 3, the piston mass (Mp) was increased by 20, 
50, 100, and 200 percent. A spring (K p ) was added to the piston 
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Figure 14. - Effect of the null center band hydraulic pump load on 
RE-1000 performance. 
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Figure 15. - Piston-displacer phase plots of the null 
center band, hydraulic-pump output Stirling engine 
for different piston masses. 
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TABLE 3. - SUMMARY OF THE DYNAMIC ANALYSIS 
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in each run sufficient to swing the additional mass at 26.7 Hz 
of case No. 5. The results are plotted in Figure 16. As the 
piston mass increased, the ratio of cyclic energy to stored 
energy became small. This means the shape of the load function 
would have less influence on the dynamics of the piston and dis- 
placer. Figure 17 clearly demonstrates this. The piston mo- 
tion became much closer to sinusoidal and the frequency 
increased. 

The dynamic analysis has been used to investigate the 
operation of the RE-1000 engine connected to a hydraulic load. 

Tne results were not intended for exact detail design, since the 
analysis itself is a simplified one. The trends, however, can 
be regarded as correct. The conclusions from the dynamic 
analysis are: 

• The system consisting of the RE-1000 FPSE and FMA de- 
signed hydraulic pump is dynamically stable. 

• At high power output, the displacer collides with the 
piston. Thin can be prevented by shortening or 
moving the piston. 

• The nonsinusoidal motion of the piston can be corrected 
by increased stored energy in the system (spring/mas3) . 

• The analysis was for the piston version of the RE-1000, 
but the trends should also apply to the flexible inter- 
face version. 


2.4.2 Thermodynamic analysis. - A number of computer simulations 
were run of the Sunpower RE- 1000 FPSE engine to determine the 
effects of increased heat transfer in the compression space re- 
sulting from the substitution of a diaphragm for the current 
piston. The simulations were done with Sunpower 's third order 
code with constrained dynamics of piston and displacer. Tables 4 
and 5 summarize the important input and output of these simula- 
tions. The code specifies a constant heat transfer coefficient 
in the compression and expansion spaces. The wall temperatures 
of the two spaces are assumed constant during the duration of an 
engine cycle. They can be held fixed or adjusted between cycles 
to approach a temperature which gives no net heat exchanae between 
gas and wall over the course of a cycle. The cylinder neat 
transfer coefficient is based on data of Taylor (ref. 1). Cylinder 
heat transfer computed in this way has given performance pre- 
dictions in close agreement with test data on *.ha current version 
of the RE-1000 FPSE. 


53 



ORIGINAL PAGE IS 
OF POOR QUALITY 



55 


gure 17. - Piston displacement versus time plot of RE-1000 with the null center 
band hydraulic load for two different piston masses. 


TABLE 4. - EFFECT OF INCREASING THE COMPRESSION SPACE HEAT TRANSFER ON THE 
PERFORMANCE OF RE-1000 FPSE WITH THE LOW PHASE ANGLE 
DISPLACER ROD 
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TABLE 5. - EFFECT OF CYLINDER HEAT TRANSFER ON THE PERFORMANCE OF 
RE-1000 FPSE WITH THE HIGH PHASE ANGLE DISPLACER ROD 
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The following is the nomenclatures for Tables 4 and 5: 


Displacer phase 
Piston swept vol 


h 

cmp 

S cmp 

T 

cmp 

h 

exp 

S 

exp 

T 


exp 

PV power 
Indicated n 


P amplitude 
cmp 

P„ mr , phase 
cmp 


Phase lead of displacer over piston 

Total volume swept in compression space 

Film coefficient in compression space 

Heat transfer surface in compression space 

Mean wall temperature in compression space 

Film coefficient in expansion space 

Heat transfer surface in expansion space 

Mean wall temperature in expansion space 

Indicated working gas power 

PV power divided by heat input to heater 

Amplitude of compression space pressure 

Phase lag of compression space pressure 
behind piston 


Regen enthalpy flux Net enthalpy transport through regenerator 

Sum pumping loss Summation of pumping power for cooler, 

regenerator and heater 


Table 4 gives results for the RE-1000 design with a 41 deg 
displacer phase angle. Table 5 gives results for a redesigned 
version of the RE-1000 FPSE with an 80 deg displacer phase 
angle. In both cases the increased heat transfer in the com- 
pression space was modeled by maintaining the identical film 
coefficient of the RE-1000 piston engine but increasing the 
surface area by the frontal area ratio of diaphragm to piston 
(9.66), or some multiple thereof, consistent with Taylor's 
correlation parameters. 

Admittedly, this approach is a crude one because it uses 
heat transfer data from ordinary internal combustion engines for 
a diaphragmed Stirling engine with quite different geometry and 
thermodynamics . 

Tables 4 and 5 show the performance of the diaphragm engine 
to be lower than that of the piston engine. This is principally 
due to a reduction in working gas pressure amplitude and phase 
lag with respect to the piston, both brought about by isothermal- 
ization in the diaphragm space. A secondary effect is the inter- 
action between the compression space and cooler due to an elevated 
average temperature in the compression space. 

A reduction in pressure amplitude is important because it 
reduces the specific power of the machine which, with nearly 
constant parasitic losses (pumping losses, regenerator enthalpy 
flu*., etc.), also reduces efficiency. The reason pressure 
amplitude decreases can be understood by noting that the pressure 
change in a gas space per piston motion is less with nearly 
isothermal conditions than with nearly adiabatic conditions. 

Calculations were made to increase pressure amplitude by 
simply increasing piston stroke (Table 4, runs 22, 23, 24). 

The pressure amplitude and power could be brought up to the 
reference level for the RE-1000 by this means but the efficiency 
was still low by a few percentage points. The reason for this 
seems to be due to the reduced phase angle of the pressure wave 
referred to the piston or diaphragm. This apparently resulted 
from a change in the temperature swing in the compression space. 
The phasing of the working gas flow through the heat exchangers 
is also affected. This can result in subtle heat transfer 
losses in the heat exchangers and reduced efficiency. It appeared 
that this problem could be overcome by optimizing the displacer 
motion in conjunction with increased piston stroke. 

Considering heat transfer in the compression space, if the 
walls of this space are considered quasi-adiabatic (Table 5, 
runs 19, 20a, 20b, 24) (no net heat exchange with the working 
gas) , the mean cycle temperature of the compression space be- 
comes higher than the cooler space wall temperature by as much 
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as 20°C. A thermodynamic loss occurs when the compression space 
gas enters the cooler. Hence, efficiency can be improved by 
active cooling of the compression space walls. Indeed, (Table 5, 
runs 21, 22, 23), when the compression space wall temperature is 
lowered to the cooler space wall temperature, efficiency rises 
somewhat . 

Even after increasing piston stroke and actively cooling 
the compression space wall, (Table 5, run 23) the efficiency of 
the diaphragm version RE-1000 could not be brought back up to 
the piston version base level efficiency (Table 5, run 20a). 

While at first surprising, this result is plausible in terms of 
the arguments presented above. 

The conclusions of the preliminary thermodynamic analysis 

are : 


• All components of the Stirling cycle must be considered 
together before any conclusions can be drawn regarding 
the merits of an increased heat transfer to approach 
isothermalization . 

• Merely increasing the heat transfer in the compression 
space of a Stirling engine is not sufficient to guarantee 
an improvement in performance. The machine must be 
re-optimized in all its aspects before near isothermal- 
ization can be expected to yield a benefit. 

• The heat transfer in the diaphragm area should be 
carefully investigated, because design and performance 
depend strongly on prediction of this phenomenon. 

• Net heat transfer from, or cooling of the diaphragm 
compression space appears necessary to achieve good 
performance . 
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3. 


DETAIL DESIGN 


3.1 Description of the Recommended Design 


3.1.1 General description. - With the preliminary design approach 
selected, it remained to produce a detail design consistent 
with the developmental program envisioned. For this the design 
incorporates : 

• Modular sections and simple, accessible parts for 
research testing. 

• Interchangeability with current RE-1000 FPSE output 
section without disturbing engine section. 

• Either vertical axis hydraulic assembly, or horizontal 
for best demonstration of dynamic balancing. 

• Displacement sensors on all engine moving parts and 
easily optional on hydraulic pump. 

• Conservative design for low stress and low cost. 

The design is described referring to Figure 18, which is a 
simplified schematic of the actual design. Basically the system 
operates as follows: 

<» Engine gas displaces oil trapped on one end of the 
piston (bottom) . 

• The piston displaces oil trapped on the other (top) 
eod, thereby displacing the bounce gas diaphragm. 

• No significant net power is transferred to this oil 
or the diaphragms, and the effective mass of this 
piston and oil is insufficient to provide the desired 
cyclic motion in combination with engine and bounce 
gas springs. 

• Additional effective piston mass is provided in the 
counter balance assembly. 

• The annular piston masses are forced by the surrounding 
trapped oil to move in opposite directions synchro- 
nously so that the entire engine system is effectively 
balanced dynamically. 
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• The annular piston mass counterbalance can be removed, 
and equivalent mass added to the piston rod, with oil 
free to flow around it, in order to demonstrate simple 
unbalanced dynamic operation similar to the normal 
RE-1000. 

• The piston rod transmits net power to the hydraulic oil 
pump, which is isolated from the diaphragm oil. 

It is obvious that the basic layout affords easy access for modi- 
fication. Detail problems in design, and the solutions offered, 
will be described piece-by-piece, from the bottom up, as it were, 
referring to detail drawings. 

3.1.2 Engine diaphragm assembly. - The engine diaphragm assembly 
interfaces with the RE-1000 mounting plate (Sunpower part: pres- 

sure vessel bottom plate) . No changes to the RE-1000 FPSE are 
required except for two through-holes in the plate for separate 
water coolant input-output, and possibly a different displacer 
rod (for improved performance, see subsection 3.2.2). 

Figure 19 is a sectional layout of the diaphragm assembly. 

A machined plug occupies the central hole in the plate vacated 
by piston version parts. This plug is the primary adapter be- 
tween the existing RE-1000 and the hydraulic output assembly, so 
it performs several functions. First, it conducts working gas 
between tho engine cooler and the diaphragm space, and it conducts 
coolant water in order to function like the engine cooler, thermo- 
dynamically. Secondly, it adapts the mechanical assembly to hold 
the displacer rod in place, to fit the existing overboard gas seal 
or O-ring, and to provide passage for the displacer LVDT sensor 
rod. Figure 20 is a two-view of this plug. The piece is made of 
aluminum, favorable for heat transfer and machining. The lower 
end has lugs to fill the otherwise dead volume of the displacer 
rod spider. The drilled gas passages are similar to those of the 
RE-1000 FPSE cooler, sized for low dead volume, flow friction, and 
cyclic thermal losses. 

Above the mounting plate is the lower half of the engine 
diaphragm support, a turned plate with milled slots. Figure 21 
is a two-view of the part. Mounting seal and shoulders suit the 
existing mounting plate. At the inner diameter is a gas seal 
with the central plug; this and the lower outer seal contain 
the coolant water. The water enters and distributes radially, 
flows semi-circumferentially in grooves around the support plate 
and central plug, and exits on the opposite radius. The plate 
is aluminum for heat transfer and machining. 

The upper surface of the plate is contoured axi-symmetrically 
to fit the extended shape of the diaphragm. In addition, there 
are radially milled small slots, tapered to zero depths outward. 
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Figure 18. - Simplified schematic of the recommended 
hydraulic output system design 
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Figure 21. - 2-View of the lower half of the engine diaphragm 
support (dimensions in United States customary units) 







to provide for gas flow distribution with low flow resistance 
when the diaphragm to support clearance is very low, that is, at 
maximum cyclic downward excursion. The "disked" contour is NC 
(numerical control) machined to specified depth versus radius 
coordinates and the unsupporting slot spans are small so that 
the diaphragm should not be damaged due to impact, sudden loss 
of engine pressure, or high oil pressure. 


The diaphragm itself is flat, annular (20.32 x 6.35 x 
0.0508 cm), and supported at its edges by clamped elastic seals, 
O-rings. This design has many advantages: 

• Structurally meridional bending and direct stresses 
are effectively nil at the edges, and the flexed shape 
offers larqe displacement for a aivei stress limit. 

• The diaphragm can be relatively thick and "springy" 
instead of "floppy" or delicate, still with moderate 
bending stresses mid span. 


• It cannot be simpler or less demanding for fabrication. 
Welding and forming are not required, only sheet stock 
having moderately high endurance limit strength (say, 
above 200 MPa for present purposes) . 


Perhaps most important, the flat diaphragm offers the 
simplest shape for control of the varying gap between the dia- 
phragm and its stroke limiting or support members. Good per- 
formance depends on achieving close clearance but free fluid 
flow; diaphragm waviness or floppiness complicates achieving 
this. Possible radial temperature gradients (less than 20°C) 
will not cause significant additional stress, and "oil canning" 
effects are far overpowered by gas-oil pressure differentials. 

The flat annular diaphragm having flexibly clamped edges, 
sometimes referred to as simply supported or articulated edges, 
has not been fully proven, but it has been used initially 
successfully in at least one Stirling engine application, by 
Cooke-Yarborough (ref. 1) , as well as for rotary seal backup 
members by FMA. Refinement of design is expected to involve 
primarily the dpt-^’ of O-ring properties and clamping groove 
section. 


The diaphragm is clamped from above, tho oil side, by an ex- 
tensively machined plate or annular ring, shown in Figure 22. In 
addition to sealing and supporting the diaphragm, it serves to 
distribute and cool the oil and, hence, t) e diaphragm, which' re- 
ceives heat from gas compression (about 7. 3W predicted, see Table 7, 
Run No. 29) . This part is also made of aluminum for good heat 
transfer and machining. The bottom surface is not only contoured 
to support the diaphragm but it also is slotted radially and grooved 
circumferentially to allow free oil flow even when the diaphragm is 
resting on it. 
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Figure 22. - 2-View of the upper half of the engine diaphragm 
support (dimensions in United States customary units) 


In addition the radial fingers formed by the slots also 
support the central plug piston of the diaphragm on the top. 

The slots, grooves, and plug piston supports are configured 
so that the diaphragm is not over-stressed even when full gas 
pressure is applied with no opposing oil pressure. 

Heat flows from the diaphragm to the fingers due to the 
cyclic convection of the oil in the immediate vicinity. Heat 
then flows through the thick fingers to the circumferential 
water cooled grooves on the upper, outer surface of the diaphragm 
support piece. The primary thermal resistance is in the oil, 
and this is practically minimal. 

The diaphragm support pieces are clamped to the RE-1000 
mounting plate by a simple steel flange ring using the existing 
bolt circle. Pressure vessel codes are satisfied with respect 
to existing bolt and seal design. 

Central to the diaphragm is a light aluminum plug piston that 
serves to support the inner edge of the diaphragm, particularly 
in the down position; to guide the cyclic gas flow; to guide the 
diaphragm in axi-symmetric displacement, without "cocking"; 
and to provide for displacement sensor rods. The diaphragm dis- 
placement is sensed directly using an LVDT rod, slug, and coil 
immersed in the oil. 

Assembly of the diaphragm parts is not unusually critical; a 
proper sequence of vertical stacking and care will avoid damage 
to the light parts such as the diaphragm, O-rings, and LVDT stems 
(e.g.: assemble diaphragm with center plug, install diaphragm 

unit and O.D. O-rings on lower support and RE-1000, install upper 
diaphragm support, upper pressure head, and last the support 
spider-ring with stationary LVDT parts) . 

The displacer motion is also sensed directly with an LVDT, 
but the rod and plug are enclosed in a thin nonmagnetic steel 
tube fixed to the plug piston and closed at the upper end. Thus 
although the displacer slug is in the gas and the LVDT coil is 
in the oil, the direct signal is not suppressed and hermetic 
separation is maintained. 

The LVDT coils are mounted on a ring having a central 
bridge that also serves to guide the plug piston rods. The 
assembly is enclosed in a simple spool piece having an electrical 
carry-through fitting. 

It is envisioned that the entire diaphragm and transducer 
assembly is assembled conveniently by sequential vertical 
stacking. 
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3.1.3 Transmission piston assembly. - Figure 23 is a drawing 
of the transmission piston and its housing, which is central 
between the engine diaphragm, bounce gas diaphragm, and balance 
piston assemblies. Operation is rather obvious: the power 
piston effectively has engine pressure on one end and nearly 
constant bounce pressure on the other, and the rod transmits 
force to the balance and pump sections. Passages are large 
relative to the piston so that flow losses and inertance of the 
transmission oil is low. Leakage past the piston and rod is 
minimized by close clearance, which is simpler and more pre- 
dictable than sealing rings for this research hardware. 

Note two items: first, the pressure force on opposite ends 

of the piston rod area itself is equalized by a flow connection 
or tube connecting the lower oil space to the top of the hydraulic 
pump assembly, as shown in Figure 18. Secondly, the output pis- 
ton assembly can be mounted with the piston rod axis horizontal, 
because the flanges for the oil flow passages are identical. With 
the piston rod horizontal, the radial facing flange of the piston 
assembly is mounted on the engine diaphragm and sensor assemblies, 
and the bounce gas diaphragm is mounted to the axial end of the 
piston assembly. Thus the bounce diaphragm appears in about the 
same place as shown in Figure 18 but t 1 •» hydraulic assembly is 
cantilevered out to the other side. '> • configuration has ad- 

vantages for demonstrating the effects, . piston dead-weight and 
inertial imbalance as compared to vertical piston motion. 

3.1.4 Bounce gas diaphragm assembly. - This assembly basically 
consists of the diaphragm, two support plates, and two flat 
flanges acting as clamps and pressure vessel, as shown in Fig- 
ure 24. The diaphragm as well as its seals and supports is 
basically the same as that of the engine. Differences are that 
there is no need for either a large central plug piston or oil 
cooling, because there is easy access to the one LVDT transducer 
and there is not substantial gas circulation or cyclic compression. 
The diaphragm should operate roughly midway between the support 
plates, and the oil side plate has circumferential grooves as 

well as radial slots to minimize flow losses. 

The gas pressure vessel attached to the assembly is optional, 
providing the correct volume. An obvious choice is the vessel 
from the RE-1000 piston version with an adapter flange. 


3.1.5 Balance {piston assembly. - The concept is not novel of 
providing inertial balance using two masses moving oppositely 
due to hydraulic coupling. The simplest arrangement for doing 
so with minimal hydraulic losses is shown in Figure 25. Both 
of the pistons, one central, one annular, have the same area; 
so their displacements are equal but opposite. Oil circulates 
only a short distance and at low velocity. The oil is effectively 
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Figure 25. — Annular piston mass counter balance assembly 
(dimensions in United States customary units) 


segregated dynamically in the upper and lower spaces due to 
small piston clearances that minimize leakage. 

Clearances and land lengths are selected to minimize viscous 
and leakage power (subsection 3.2) but there will be enough 
leakage past the pistons and rod bushings so that the annular 
piston will drift in either direction during operation, and down- 
ward after shutdown. An oil bleed compensation system is nec- 
essary to locate and stabilize the positions of all the pistons 
and diaphragms, to be discussed in subsection 3.1.7. 

Note that the balance piston is not necessary for engine 
operation. Initial testing of the engine and vertical pump rod 
configuration should be carried out using an appropriate single 
small-diameter mass clamped to the rod in the balance section. 


3.1.6 Hydraulic pump system. - This system consists of the rod 
having a small integral piston, the pump block, two sets of in- 
put and output check valves, small accumulators local to the 
block, a pumped storage accumulator, and an oil sump. All but 
the last two items, which can be remote and sized optionally, 
are shown in Figure 26, with several sectional views. The design 
and operating principles of the pump are basically the same 
as almost all double acting piston pumps having check valves, 
except for the null center band. As described previously 
this is the simplest, most inherently stable way to start and 
run a FPSE coupled to an efficient hydraulic output. 

However, this is not without risk of penalty, depending on 
design. The risk is caused by the fact that at the beginning 
or the end of the null band, the piston has significant velocity, 
causing a sudden demand for intake and discharge flows to stop 
or start, respectively. (In a normal pump having harmonic pis- 
ton motion, valve action is mostly at low piston velocity, TDC 
and BDC , and flows accelerate smoothly) . Consequently, the 
present pump is designed with large flow passages, light valves, 
and minimum passage lengths to the attached accumulators, all 
of which tend to accommodate the "spurts" of flow with minimal 
hammer and cavitation. 

The pump is constructed of a monolithic block having drilled 
passages with inserted plugs and valve modules to direct the 
flow. Bushings trapped by the end caps guide the piston and rod, 
which are one piece. Accumulators and valve modules are acces- 
sible without dismantling the pump. The discharge valve is a 
typical ball-check, sprung and guided. The intake valve is a 
"thumbtack" shape to achieve the necessary orifice area and yet 
lightness to admit flow from an atmospheric sump without de- 
structive cavitation. 
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Figure 26. - Hydraulic pump system (dimensions in United States 
customary units) 



The cross ports are drilled such that the piston is always 
guided and laterally pressure-balanced in its bore, and there 
is ample overstroke distance. In addition, the drilled ports 
for null band, internal-bypass oil circulation provide a some- 
what gradual cutoff or opening as the piston occludes a diminish- 
ing port sector. Otherwise, leakage past the piston is minimized 
by close clearance, without sealing rings. 

The pump as shown is sized for 2kW output with mechanical 
and volumetric efficiencies of 80 percent and 14 MPa pressure 
differential. This is the maximum power expected of the engine. 
If less power is available at the expected stroke, it can be 
accommodated at reduced pressure differential. Otherwise, the 
engine balance point will be reached at shorter stroke and even 
less power. A prudent plan may be to start with a smaller pis- 
ton and cylinder bore and rework the parts as the engine is 
optimized . 


3.1.7 Control methods . - Control of this, or probably any 
hydraulic output FPSE, involves two general aspects: startup 

and stabilization of power; and location and stabilization of 
all three or four displacement components, namely the two dia- 
phragms and one or two piston balance masses. The system has 
at least two more degrees of freedom than the "simple" FPSE, 
so it requires more design effort and system means to ensure 
stable operation at the design parameters. We reiterate the 
obvious: many design parameters are not "cast in iron" with a 

FPSE, and these desired parameters will not be obtained or remain 
constant unless they are forced to. 

The simplest aspect is startup and control of power, 
yiven proper average fluid distribution parameters such as 
pressures, volumes and temperature. Then, when stopped, both 
diaphragms and the power piston will be near mid-stroke and the 
displacer will rest on bottom, in the hot end. For startup, 
using suitable switching valves, first the oil pressure is bled 
from the hydraulic piston rod end, raising the piston and 
positioning the engine diaphragm to begin a compression stroke. 
Figure 27 is a schematic of the engine system (from Figure 18) 
showing the necessary valve functions. When pressure is relieved 
from the upper rod end, engine gas-oil pressure will cause the 
piston to rise, because of net-force imbalance. After any short 
dwell time, the displacer should be at hot or bottom position. 

Then if the bleed valv' . are switched back to normal, then the 
engine should start oi * ae rebound downward compression stroke 
of the piston. Better, if the bleed is shut and higher pressure 
accumulator oil is introduced to the rod end then a more forcible 
rapid compression stroke will occur. In this case, the high pres- 
sure must be switched off and the rod end oil bypass opened quickly 
before the end of the compression stroke. It appears likely that 
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either of these techniques will result in proper phasing of the 
displacer, and there is no load on the engine within the null 
band of the pump, so "hot" startup should be easy. 

"Cold" starts should be accomplished in the same way, except 
that all moving parts must be brought to proper initial positions, 
and the hot end of the engine raised to a sufficient temperature. 
The most likely state of the engine after an extended shutdown is 
with the diaphragms resting on their oil-side supports and no oil 
pressure. A bleed-fill system is fairly obvious that it can be 
used during research testing, operator controlled with position 
references, the LVDTs, in order to restore oil pressure and place 
the two diaphragms in good positions intermediate to their sup- 
ports. Both the output piston and the counter-balance piston will 
have sunk to bottom positions before a cold start is attempted. 

The piston still can b • lifted by bleeding oil from the upper rod 
end, but the motion will be slow, controlled by leakage through 
the counter-balance annuli. A more automatic cold start sequence 
may be conceived after a system to maintain mean-position sta- 
bility is proven; this is to be discussed later. 

After the engine is running, given stable, average or 
center positions of the moving components, power control is 
straightforward: the existing system controls heater temperature 

and the hydraulic stroke grows until the available power is 
absorbed, as described before for the null center band pump con- 
cepts. Change in power level is not normally desired while 
filling the output hydraulic accumulator, but it can be provided 
for test purposes by changing the engine temperature, the charge 
pressure, the mean dead volume (by diaphragm mean position), or 
the hydraulic accumulator pressure, leading to different strokes 
at different output levels. Shutdown is simply a matter of 
turning off the heat, and obviously this can be signalled by an 
oil pressure or level sensor. 

As implied earlier, keeping the diaphragms and pistons from 
drifting away from their mid positions is the "tricky" part. 

Not counting the engine displacer, the system has three degrees 
of freedom with the two diaphragms and the piston assembly. The in- 
ertial counter balance adds a fourth degree of freedom if it is 
used. In this system there are two forcing functions independ- 
ently related to separate displacements in the engine and the 
hydraulic pump. For comparison, the piston version RE-1000 has 
only one degree of freedom, the piston (again not counting the 
displacer) , and only one independent forcing function, which is 
related to piston displacement. The piston acts directly on the 
engine working gas as well as the gas in the output damper. In 
the hydraulic output version the engine gas displacement is sep- 
arated from the hydraulic pump displacement by an independently 
varying (i.e.; unknown inter-chamber leakage) volume of oil that 
acts to transmit force and motion. 
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The only possible control advantage of the hydraulic dia- 
phragm version is that there is no leakage between the gas spaces 
except that provided in a control system. On the other hand, 
there is uncertainty concerning many factors: the directional 

bias forces of the diaphragms as they approach their supports, 
the bias force leading to drift of the hydraulic piston, and 
the net leakage of the pump oil to the transmission oil. With 
respect to first order effects, there are neither bias nor centering 
forces in the output transmission system itself, except for the 
piston weight, and there is no oil leakage bias (7 MPa transmission 
oil, 14 MPa accumulator, open sump) . The primary bias, the piston 
weight, can be balanced permanently with the appropriate additional 
pressure (charge mass) in the engine versus that in the bounce 
space . 

Consequently, the tentative conclusion is that with the 
engine and bounce gas hermetically isolated, the positional 
drift of all parts will be slow enough for operator response. 
Referring to the three position sensors, we recommend the follow- 
ing control method for initial testing: adjust the gas balance 

in order to control engine diaphragm position and balance the 
vertical piston force; adjust the balance of oil on either side 
of the transmission piston in order to center the output piston; 
inject or bleed oil from the balance space in order to maintain 
transmission oil volume and, hence, bounce gas diaphragm position, 
given that the engine diaphragm mean position is constant. Engine 
diaphragm position will /ary primarily only with engine tempera- 
ture, or effective "reduced" volume per mass. A system of bleed- 
check valves to provide these control functions is shown in 
Figure 28. This system has not been sized because the prediction 
of drift due to obscure secondary effects is difficult. 

Systems to automatically control the positions of the hydrau- 
lic transmission parts have been conceived but not designed in 
detail. The tentatively preferred scheme is shown in Figure 29. 

It consists first of oil bleed valve spools attached to both 
diaphragms; they dump oil when a diaphragm approaches too near 
its gas-side support, thereby preventing impact and providing 
continuous operation with minimal dead volume in the engine gas 
compression diaphragm space. Secondly, a limited amount of oil 
flows from the high pressure accumulator into the cylinder annulus 
around the output piston. Bleed passages between the oil on 
opposite sides of the piston are modulated by the piston position 
so that a time-average bias in the leakage flow tends to center 
the piston. (This "centering-port" method is often used with gas 
leakage to stabilize free piston engines internally.) Similar 
bleed-bias ports are required to maintain approximate centering 
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Figure 29. - Schematic of automatic position control 
for hydraulic output components 
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of the balance-piston motion in its cylinder. Third, it appears 
necessary to balance the engine gas and bounce gas pressures in 
order to avoid eventual impact or "bottoming out" of diaphragm 
or piston. This can be done by using the oil flow dumped from 
the diaphragm oil valves to adjust differential gas bleed between 
engine and bounce spaces. For example, a gas spool valve actuated 
by oil pressure through small di phragms could be incorporated. 

Other alternatives may be considered for the automatic system. 
For example, a spool valve attached to the gas side of the 
bounce diaphragm could replace the oil dump activated gas bleed 
valve components. Simply but mechanically, a stiff spring sus- 
pension of the transmission piston might eliminate the need for 
any gas bleed system if a suitable gas charge distribution could 
be set and held. 

Finally, if a counter-balance piston is used, two bleed 
ports cyclically uncovered by the counter-balance annular piston 
will provide centering by increasing the average pressure in the 
lower half of the counterbalance oil filled chamber. Because 
this pressure will also exactly balance the transmission piston 
weight, this arrangement will be insensitive to orientation. 
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3.2 Detailed Analysis 


This section describes the analyses performed in parallel 
with detail design of the hydraulic output system for the RE-1000 
FPSE. The detailed analyses include the metal diaphragm analysis, 
the thermodynamic and dynamic analyses of the engine, and esti- 
mates of additional power losses due to hydraulic output con- 
version of the RE-1000 FPSE. 


3.2.1 Metal diaphragm analysis. - Two metal diaphragms are used 
in the hydraulic output version of RE-1000 FPSE. They are the 
engine diaphragm and the bounce gas diaphragm. The engine dia- 
phragm provides a hermetically sealed interface between the engine 
gas and the hydraulic fluid. The bounce gas diaphragm separates 
the hydraulic fluid from the bounce space gas. 

The diaphragm analysis will provide the information on the 
following items of interest: 

• Stress level and distribution - this information is 
required for the design of the edge sealing and sup- 
porting method. It will also be used to ensure that 
the maximum stress of the diaphragm remains well below 
the endurance limit of the diaphragm material. 

• Axial, displacement distribution - this information will 
be used to design the over-pressure protection sup- 
ports. The squeezing and trapping of gas or oil be- 
tween the diaphragm and the supports should be pre- 
vented by careful contouring of the support structures. 

• Meridional displacements of edges - this is critical 
information for design of the edge sealing and sup- 
porting. Especially, if the O-ring seals are used for 
the edge sealing, the meridional displacement of edges 
will cause the O-rings to roll back and forth radially. 
This rolling distance should be short to prevent the 
fatigue failure of the O-ring material. 

• Stiffness of diaphragm near TDC and BDC - the program 
NONLIN can analyze only static loading case. The 
dynamic effect was approximately accounted for by 
introducing a static force applied at the inner edge. 

The magnitude of the force was set equal to the inertia 
force of the center plug mass oscillating at design 
point conditions. The stiffness of the diaphragm (k) 
and the center plug mass (m) introduce a natural 
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frequency equal to ( 1 / 2tt ) / k/m . In order to prevent 

(high frequency oscillation) of the diaphragm at its 
TDC and BDC , the above natural frequency should be 
designed to be much higher than the engine frequency. 

For the analysis of the diaphragms, two predeveloped com- 
puter programs, NONLIN and ANSYS, were utilized. The program 
"NONLIN" is a bellows/diaphragm analysis developed by Dr. J.F. 
Lestingi, et al. (ref. 2) at Battelle Memorial Institute, Columbus 
Ohio. It is based on a multisegment numerical integration of a 
set of simultaneous nonlinear first order differential equations. 
This program is capable of analyzing linear elastic axisymmetric 
and nonsymmetric deformations, and nonlinear axisymmetric de- 
formations of bellows and diaphragms. 

The "ANSYS" program is a proprietary, general purpose finite 
element program developed by Swanson Analysis, Inc. The program 
is available through a commercial computer network. 

FMA obtained the program NONLIN through the courtesy of 
Battelle Institute, Ohio. NONLIN was used for the actual analysis 
of the diaphragms. ANSYS was used only to check the results 
from NONLIN in the initial stage of the analysis effort. 

Initial sizing of the diaphragms was done by approximate 
calculations combining the Belleville spring analysis (ref. 3) 
and simple bending theory for flat plates. Table 6 shows the 
specifications and boundary conditions of the two diaphragms 
used in the analysis. The specifications given here represent 
the final diaphragm designs for the hydraulic output RE-1000 FPSE. 

The results of the diaphragm analysis are presented graph- 
ically in Figures 30 through 37. Brief discussions of each follow 

Figure 30 shows the diaphragm center deflection as a function 
of the pressure loading on the diaphragms. A flat annular dia- 
phragm has a very low stiffness near the mid position. The stiff- 
ness increases rapidly as the diaphragm is deflected away from the 
mid position. Beyond a certain point, the diaphragm stiffness be- 
comes almost constant as shown in Figure 30. Because of the 
rapidly increasing stiffness, the program NONLIN required a very 
small load step in order to achieve convergence. A 0.254 mm 
thick diaphragm was initially analyzed. As shown in Figure 30, 
the 0.254 mm diaphragm had such a low stiffness that only 1380 Pa 
pressure differential was required for the design point displace- 
ment of 80 cm . The stiffness in this case was so low that even 
a small radial pressure gradient due to fluid squeezing would 
cause the diaphragm to deform locally in a totally unpredictable 
way. That is why the thickness was increased to 0.508 mm in which 
case the diaphragms were stiff enough to avoid unwanted local 
deformations . 
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TABLE 6. - SPECIFICATIONS AND BOUNDARY CONDITIONS 
FOR THE DIAPHRAGM ANALYSIS, NONLIN 


Specifications 

Engine diaphragm 

Bounce gas diaphragm 

Outer diameter (cm) 

20.32 

20.32 

Inner diameter (cm) 

6.35 

2.54 

Thickness (cm) 

0.0508* 

0.0508 

Young's modulus (MPa) 

2.069 x 10 5 

2.069 x 10 5 

Poisson's ratio, V 

0.27 

0.27 

Boundary conditions 



Outer edge 

Simply supported to the 
casing 

Simply supported to the 
casing 

! Inner edge 

Simply supported to the 

Axially clamped to and 

| 

center plug which moves 
with the diaphragm 

radially guided by the 
center plug which moves 
with the diaphragm 

Volume amplitude at 



design point (cm3) 

40 

40 
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Figure 32. - Circumferential stress distribution of the engine 
diaphragm at design point. 
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gure 33. - Meridional stress distribution of the engine 
diaphragm at design point. 
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Figure 35. - Meridional stress distribution of the bounce gas 
diaphragm at design point. 


OUTER SURFACE 






The spring constants, k, for both diaphragms near the TDC 
or BDC positions were about 96,000 N/m. For the engine diaphragm 
with the center plug mass, m, of 181g, the natural frequency, 
f, obtained from the equation 



is close to 120 Hz. For the balance diaphragm with the center 
plug mass of 45g the natural frequency is close to 240 Hz. 

In each case, the natural frequency of the center mass sprung 
by the diaphragm spring is much higher than diaphragm oscillating 
frequency of 30 Hz. Therefore, the high frequency ringing of 
the diaphragm near TDC or BDC will have very small amplitudes 
so that it will have very little effect on the overall diaphragm 
dynamics . 

Figure 31 shows how much the engine diaphragm moves axially 
and meridionally under a static pressure load which gives the 
design volume displacement. The axial displacement curve is 
used to get the right shape for the supports. The volume dis- 
placement was numerically calculated from the axial displacement 
curve. The meridional displacement curve indicates that quali- 
tatively speaking, the inner half of the annular diaphragm moves 
radially outward, and the outer half moves radially inward. The 
distances the edges move radially are very small (< 0.03048 mm) 
which will cause little rolling fatigue problem for the O-ring 
seals at the edges. 

Figure 32 gives the circumferential stress distribution in 
the engine diaphragm at design point. It shows that the maximum 
circumferential stress of 103 MPa occurs where the diaphragm 
is connected to the center plug. About two- thirds of the maxi- 
mum stress is due to the membrane stretching and about one-third 
is due to the bending. 

Figure 33 shows the meridional stress distribution in the 
engine diaphragm. It clearly shows the advantages of simply 
supported edges: the meridional stresses are zero at the edges. 

Figures 34 to 36 snow the distribution of displacement, cir- 
cumferential stress, and meridional stress in the balance dia- 
phragm. As shown in Figure 34, the distance that the outer edge 
moves is less than 0.03048 mm as was the case with the engine 
diaphragm. Since the inner edge of the diaphragm is clamped 
axially to the plug, the meridional stresses in Figure 35 show 
considerable level of stresses, up to 48 MPa at the inner edge. 
Still the maximum stress in the balance diaphragm is only about 
90 MPa in the circumferential direction as shown in Figure 36. 
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Figure 37 shows the stress versus cycles to failure data for 
AM 350 SS sheet of 1.016 mm thickness. Although our diaphragm 
design has half the thickness of the fatigue test specimen, the 
endurance limit is not expected to depend significantly on thick- 
ness within practical design ranges. Tho maximum equivalent 
stresses of the two diaphragms are 100 to 110 MPa which is well 
below the endurance limit of 590 MPa for AM 350 SS. 


3.2.2 Thermodynamic Performance Predictions 


3. 2. 2.1 Methods of analysis. - For the detailed analysis of the 
hydraulic output Free Displacer Stirling engine, Sunpower used 
their third order Stirling engine analysis. 

There are two modes in the above analysis: the constrained 

mode and the unconstrained mode. The constrained mode deals 
only with the thermodynamic aspect of the Stirling engine. The 
dynamics of the reciprocating components are predetermined by 
a separate dynamics analysis using a simplified thermodynamics. 

The unconstrained mode combines the dynamics of the piston, the 
displicer and the load, and thermodynamics. This mode takes into 
account the interaction between the dynamics of the moving parts 
and the engine thermodynamics by solving them simultaneously. 

The unconstrained mode requires many cycles of calculation in 
order to reach the steady state, and therefore is quite expensive 
to use. 

In the present problem, the null center band hydraulic pump has 
a load function which resembles truncated square waves. The highly 
impulsive or step function like nature of the load necessitated 
the use of the unconstrained mode to closely simulate the actual 
processes of the engine. The less expensive, constrained mode 
analysis was used to obtain the performance trend and to perform 
a sensitivity analysis of important operating parameters. 


3. 2. 2. 2 Critical parameters. - Within the scope of the design 
program, there are two parameters significantly affecting the 
power output and the efficiency of the engine. They are displacer 
phase angle with respect to the piston motion and the cylinder 
heat transfer. Brief discussion of each follows. 


3. 2. 2. 2.1 Displacer phase angle. - The displacer phase angle 
is closely linked to the overall engine dynamics. Varying the 
displacer phase angle will change the frequency, the amplitudes 
of the piston and the displacer, the pressure wave phase angle 
and also the engine performance. It will also change the fre- 
quency slightly. For a displacer of a given mass, the phase 
angle depends on the following: 
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• The displacer gas spring constant which is inversely 
proportional to the mean volume of the gas spring. 

• The displacer load which is determined by the pressure 
drop loss through the heat exchangers and the gas spring 
hysteresis loss. 

The phase angle of the displacer can be easily changed by 
varying gas spring volume. The original RE-1000 FISE had two 
displacer rod sets, the main difference of which was the gas 
spring volume. The rod set with the larger gas spring volume 
had the low displacer phase angle of 45 deg, which was determined 
by the optimization of the engine for the highest efficiency. 

The rod set with the smaller gas spring volume had the displacer 
phase angle of 80 deg. The high phase angle would result in th® 
highest power. 

For the hydraulic output Stirling engine, however, the un- 
constrained mode analysis indicated that the rod set should be 
redesigned as will be discussed later. 


3. 2. 2. 2. 2 Cylinder heat transfer. - With the increased heat 
transfer area and the extremely narrow spaces in the compression 
space of the diaphragm engine, the magnitude of the heat trans- 
fer coefficient has a significant consequence in the predicted 
performance of the engine. 

Two methods were used to estimate the heat transfer c fficient 
inside the cylinders: Taylor's Internal Combustion (I.C.) engine 

correlation (ref. 5) and FMA's cyclic heat transfer analysis 
(ref. 6, 7). 

Taylor's I.C. engine correlation is based on data gathered 
from 19 different types of engines for automobiles, aircrafts, 
trucks, locomotives and other commercial use. It is expressed 
by 


k .0.75 

h = 10.4 -2 ' GB 

D 


( — ) 
\W 0 / 


( 2 ) 


where 


h = 
k = 

g 

B = 


Heat transfer coefficient, based on piston top area 

Conductivity of the engine gas 

Bore 
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G = Average mass flow rate per cycle divided by piston area 

y = Viscosity 

g » Gravitational constant. 

Sunpower previously observed that when the above correlation 
had been used in their third order analysis, it had resulted 
in a close agreement of the analysis with the test data. Also, 
for normal bore/stroke cases, the values of the heat transfer 
calculated from the above equation were not too different from 
the values calculated by FMA's cyclic heat transfer analysis 
which is described in the following paragraphs. 

It is fair to say that the Taylor's formula is a very 
valuable empirical correlation for the heat transfer in cylinders. 
But the lack of combustion process and values in tho Stirling 
engine makes it difficult to be accepted without actual testing 
of Stirling engines. Especially, this correlation is invalid 
for high bore/stroke case such as diapnragm piston engines. 

The FMA's cyclic heat transfer analysis has been partially 
published (ref. 6, 7). The analysis is valid for 

• The normal bore /stroke cylinder-piston cases. 

• The high bore/stroke cylinder-piston or close gap 
diaphragm cases. 

For the normal bore/stroke case, the FMA analysis utilizes 
empirical heat transfer enhancement factor correlation based on 
an unpublished gas spring data (ref. 6). 

For the narrow gaD case, the FMA analysis includes the 
effect of mixing or eddy conductivity. Although no test data 
exist, the effect of mixing modeled here reoresents a conser- 
vative estimate due to the fact that a steady-flow, not periodic 
flow, assumptions were used to obtain the heat transfer enhance- 
ment factor. 


3.2.2. 3 Results of the third order analysis. - The results 
include: 

• The base-point, or reference simulations of the stand- 
ard RE-1000 FPSE with the 80 deg phase displacer and 
the diaphragm piston engine. 

• The sensitivity analysis of the diaphragm piston 
Stirling engine with respect to the cylinder heat 
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transfer h, the phase ?.ngle of the displacer, $<5# and 
the pump load force, F- oac j- 

• The predicted performance of the diaphragm piston 
Stirling engine. 


3.2.2. 3.1 Reference simulations. - Sunpower used the constrained 
mode third order analysis for the reference simulations. The 
purpose of the reference simulations was to establish the base- 
point performance of the standard RE-1000 FPSE to be compared 
with the predicted performance of the diaphragm piston engine. 
Also, the constrained mode simulation of the diaphragm piston 
engine serves as an initial condition for the unconstrained 
mode simulation. 

In each simulation, FMA supplied the heat transfer coeffi- 
cients for both expansion and compression spaces from the cyclic 
heat transfer analysis mentioned earlier. 

Run 20a and run 28 of Table 7 show the results of the simu- 
lation of the standard RE-1000 FPSE with Taylor's heat transfer 
correlation and FMA heat transfer correlation, respectively. 

Run 29 gives the results for the diaphragm engine simulation. 

Note that the piston swept volume for the diaphragm engine is 
36 percent higher than that of the standard RE-1000 FPSE. The 
increased stroke was to compensate for the decrease in mass 
flow amplitude for the near isothermal diaphragm engine, and to 
maintain the pressure amplitude level. 

The displacer rod is of the high phase angle design (80 deg) 
for all three simulations. 

The magnitudes of the h eat transfer calculated by Taylor 
correlation are 14 to 118 percent higher than those predicted 
by FMA correlation. This can be explained by the fact that the 
I.C. engines have a violent combustion period and an intentionally 
augmented inlet turbulence level which is known to persist over 
the entire cycle whereas the Stirling engines do not experience 
either condition. 

Although the increased piston swept volume maintains the 
pressure wave amplitude, the increased heat transfer in the com- 
pression space of the diaphragm engine changes the engine per- 
formance noticeably. It considerably lowers the pressure wave 
phase angle and more than doubles the available energy loss in 
the compression space. The end result is a drop in power and 
efficiency as shown in Table 7. In run 29, the compression space. 
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TABLE 7. - RE-1000 REFERENCE SIMULATIONS CONSTRAINED THIRD ORDER ANALYSIS 
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especially the diaphragm area experienced a large amount of heat 
transfer. Under the assumption of the constant wall temperature 
at the sink temperature, the compression space had a net heat 
transfer of 1370W out of the space. About half of this heat 
had to be transferred through the hydraulic side, which neces- 
sitated the heat transfer augmentation measures in the hydraulic 
side of the diaphragm described in subsection 3.1. 


3. 2. 2. 3. 2 Sensitivity analysis. - 

Heat transfer in the compression space: From run 29 of 

Table 7 the large available energy loss in the compression space 
of the diaphragm engine reflects the fact that the processes are 
far from being isothermal. Here, the available energy loss refers 
to the decrease in available portion of tne heat due to heat 
transfer across finite temperature differences and due to mixing 
of gases at different temperatures. It would be interesting to 
find out how much heat transfer increase or decrease would cause 
an appreciable improvement in performance. 

Figure 38 illustrates the effect of varying the magnitude 
of the heat transfer in the compression space on the predicted 
performance of the diaphragm engine. Included in Figure 38 
are PV work, indicated efficiency and available energy loss for 
the compression space as a function of the heat transfer (hA) in 
the compression space. The analysis was done using the Sunpower 
constrained mode third order analysis. The detailed operating 
conditions and performance parameters are listed in Table 7. 

From Figure 38 the following observations were made: 

• The P-V power of the adiabatic end (small hA) is higher 
than the isothermal end (large hA) . This is because 
the mass flow in and out of the compression space for 
the near adiabatic case is higher than that of the 
near isothermal case. The increased mass flow ampli- 
tude is translated into an increased pressure wave 
amplitude, and thus a higher power output as shown in 
Table 8. 

• The efficiency for the near isothermal case is much 
higher than that of the near adiabatic case. It is 
explained by the difference in the available energy 
loss in the compression space for the two cases. In 
the perfectly isothermal case, the available energy 
loss would be zero. In the perfectly adiabatic case, 
the available energy loss would not be zero because 
there is still the mixing of gas at different temper- 
atures in the compression space* 


® RUN 32. RE-1000, DIAPHRAGM, NEW RECONMENDED 
ROD, UNCONSTRAINED DYNAMICS INCREASED 
COMPRESSION SPACE DISPLACEMENT 
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Figure 38. - Predicted performance of the diaphragm piston 
free displacer Stirling engine as a function of heat 
transfer in the compression space. 
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TABLE 8. - THE PERFORMANCE OF HYDRAULIC OUTPUT, FREE 
DISPLACER, DIAPHRAGM PISTON STIRLING 
ENGINE VERSUS COMPRESSION SPACE 
HEAT TRANSFER 


— Run No. 

Parameters* 

26 

27 

29 

Displacer phase (deg) 

80.2 

80.2 

80.2 

Piston swept volume (m^) 

80.1E-6 

80.1E-6 

80.1E-6 

h (W/nT-K) 

cmp 

1. 36E5 

4.5E3 

3.07E4 

h S (K/K) 

cmp 

1.16E4 

3.83E2 

2.61E3 

T (K) 

cmp 

313 

313 

313 

h (W/m 2 -K) 

exp 

4.88E3 

4.88E3 

4.88E3 

h S (W/K) 

exp 

6.84E1 

6.84E1 

6.84E1 

T (K) 

exp 

835 

822 

827 

PV power (W) 

2210 

2700 

1930 

n i. 

0. 384 

0.366 

0.311 

P amplitude (N/m 2 ) 

cmp 

1.25E6 

1.65E6 

1.35E6 

P phase (deg) 

cmp 

21.4 

18.6 

16.0 

Regenerator enthalpy 
flux (W) 

1080 

1060 

1020 

Sum pumping loss (W) 

364 

409 

338 

Available energy 
loss in the compression 
space (W) 

210 

570 

740 


*For the nomenclature, refer to subsection 2.4. 
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• In the intermediate reqion between the adiabatic and 
the isothermal case, the efficiency and the power drop 
to their respective minimum values. This is because 
the available enerqy loss peaks in the intermediate 
reqion as shown in Fiqure 38. As it happens, the run 
29 of Table 7 is in this fairly wide, low performance 
reqion. An interesting thinq about this is that even 
if our estimate of the heat transfer was off by a 
factor ol 2 in either direction, the performance of 
the diaphraqm enqine 'would not differ much. 

• Toward the adiabatic end, the piston enqine (run 28,) 
puts out much loss P-V power than the diaphraqm enqine 
(solid curve), although the efficiencies seem to match 
approximately. This is only because the diaphraqm 
enqine had 36 percent hiqher piston swept volume than 
the piston enqine, and beinq adiabatic, there is no 
decrease in effective stroke as with a near isothermal 
case . 

Dis placer phase angle : Using the constrained third order 

analysis, the performance of the diaphraqm engine was simulated 
tor a wide range of the displacer phase angle. Table 9 gives 
the detailed information on operating conditions and output 
variables. The Taylor's correlation (ref. 5) was used for the 
heat transfer calculation. The increased piston swept volume 
was used as before. The pressure wave amplitude increases as 
the displacer phase anqle increased from 40 deg to 100 deq. 

The phase angle* of the pressure wave has a maximum near 70 deg. 

The P-V power, which is directly related to the amplitude and 
the phase of the pressure wave peaks around 80 deg of displacer 
phase angle. Figure 39 shows the efficiency and power variation 
with respect to the displacer phase angle. The maximum efficiency 
occurs around 60 deg. Note the shift in maximum efficiency 
phase angle: the original RE-1000 FPFE had been designed to 

have 45 deg phase angie for maximum efficiency. The 15 deg phase 
shift is thought to be due to the increased heat transfer in the 
compression space of the diaphragm engine. 


Hy draulic p um^> load force: The dynamically unconstrained, 

third order analysis was used to predict what would happen if 
the hydraulic load force was varied over a range of 900 to 1300 
Newtons. The following is the brief description of how the null 
center band hydraulic pump is incorporated in the unconstrained 
third order analysis. 
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TABLE 9. - THE PERFORMANCE OF DIAPHRAGM PISTON, FREE 
DISPLACER, HYDRAULIC OUTPUT STIRLING 
ENGINE VERSUS DISPLACER PHASE ANGLE 


Run No. 

25-1 

25-2 

25-3 

25-4 

25-5 

25-6 

25-7 

Displacer phase (deg) 

40 

50 

60 

70 

80 

90 

100 

Piston swept volume (m 3 ) 

80. 

IE-6 

h (W/m 2 -K) 

cmp 

- 

1 . 

36E5 (Taylor's 

correlation) 

S ^ (m 2 ) 
cmp 

2. 

48E-2 


T (K) 

cmp 

313 

h (W/m 2 - K) 

exp 

5. 

68E4 


S (m 2 ) 

exp 

2. 

57E-3 



T (K) 

exp 

823 

822 

822 

825 

830 

835 

843 

PV power (W) 

1185 

1534 

1810 

1993 

2056 

2036 

1947 

Indicated 

efficiency 

0.315 

0.336 

0.344 

0.341 

0.326 

0.310 

0.293 

AP (N/m 2 ) 

cmp 

9.54E5 

1.05E6 

1.15E6 

1. 26E6 

1 . 37E6 

1 . 46E6 

1 . 56E6 

Ap c*p phase 
(deg) 

a 

16.7 

17.9 

18.0 

17.1 

16.0 

14.5 

Regenerator . 
enthalpy flux (W) 

673 

745 

850 

945 

1049 

1130 

1270 

Sum pumping 
loss (W) 

168 

204 

249 

302 

364 

427 

497 















































































ORIGINAL P<»-- 
OF POOR QUALITY 


The hydraulic load is described mathematically by a sub- 
routine which offers a square wave force, F to the engine piston 
defined as follows: 

F = 0.0 if abs (x ) < 1 . 1 cm 

p (3) 

F “ ± F LOAD if abs (X P } l Llcni 

where the sign of Fload depends on the piston velocity, and 
Xp refers to the displacement amplitude of the hydraulic trans- 
former piston. In other words, there is no force for a null 
center band around the piston mid-position of 55 percent of the 
design stroke of 4.0 cm. Outside the null center band, the 
piston sees a force of + FlqaD' w ^ ere the positive force opposes 
the piston motion. Fload the P r oduct of hydraulic pressure 

and the pump piston area and is chosen to absorb the correct 
power from the engine at the design point. 

The diaphragm engine simulated for this purpose has the 
high phase angle displacer rod set, the effective piston mass 
of 2.8 kg, and a fixed null center band of 2.2 cm. Figure 40 
shows that as the load force increases, the load clearly shows 
the tendency of restricting the piston motion: the piston 

stroke, the frequency, and the power decrease by 30 percent 
and 43 percent, respectively as the load force is increased 
from 900 to 1300 newtons. The displacer phase angle remained 
almost constant despite all the other changes. In these sim- 
ulations, there was no transient overstroking or collisions 
which again demonstrates the inherent dynamic stability of the 
oresent design. 


3. 2. 2. 3. 3 Performance prediction for the diaphragm-piston, 
f ree-displacer , hydraulic output Stirling engine. - The diaphragm 
engine was to be designed to have the displacer phase angle of 
80 deg for maximum power production. Therefore the reference 
simulation run 29 of Table 7 shows 80 deg displacer phase. How- 
ever a dynamic analysis performed on the diaphragm enaine cf 
run 29 indicated that the gas spring stiffness in the existing 
displacer spring of the RE-1000 FPSE was too stiff to give de- 
sired displacer phasing in the unconstrained mode. The dynamic 
analysis suggested a spring stiffness of 2.95E4 N/m whereas the 
existing spring has a stiffness of 4.15E4 N/m. It was decided 
to simulate the existing gas spring. A s^ies of simulations 
were made with various piston masses and hydraulic discharge 
pressures to obtain an acceptable run. 
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UNCONSTRAINED DYNAMICS 
HIGH PHASE DISPLACER ROD 
EFFECTIVE PISTON MASS 2.8 Kg 
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The simulations were started with initial conditions pro- 
vided by the output of constrained run 29. In other words the 
piston and displacer were already moving with design point 
amplitude and phasing. All simulations were dynamically stable, 
either decaying or increasing slightly in amplitude to a steady 
operating cycle. No evidence was seen of runaway amplitude over 
the 40 percent variation in hydraulic discharge pressure 
investigated . 

Simulations 32d and 32e described in Table 10 gave the best 
dynamic behavior, however, efficiency for both was low due to 
improper displacer phasing as a result of. the stiff displacer 
spring . 

Figure 41 shows a piston-displacer phase plot of ru.i 32d. 

The interaction of the engine with the non-linear hydraulic 
pump load produces a slightly egg-shaped phase diagram instead 
of a perfect ellipse obtained for sinusoidal piston and displacer 
motion . 


S uggested displacer spring modificatio n: Since modification 

of the displacer gas spring stiffness is a relatively simple 
shop operation, it was decided to simulate the system with the 
correct gas spring stiffness of 2.95E4 N/m. 

A sketch of the existing displacer rod and a modification 
which would increase the gas spring volume and hence reduce 
its stiffness is shown in Figure 42. The modification increases 
the gas spring volume from 17.9 cm3 to 25.1 cm3 by boring out 
the rod and inserting a cup liner. The centering leak flows 
through the annular gap between the cup and the rod. The pres- 
sure pick up channel is isolated from the centering leak by a 
tube segment which connects the interior of the gas spring to 
the existing hole in the displacer rod. 

Run 33 of Table 1C summarizes the results of this run. 

The efficiency is good as is the power level. The piston stroke 
is a bit too high which could be brought down by increasing 
the hydraulic pump discharge pressure slightly. 

Table 11 shows the combined results of the engine thermo- 
dynamic performance predictions given in Table 10, and the 
additional loss estimates of subsection 3.2.3. Table 11 clearly 
recommends the use of the new revised displacer rod shown in 
Figure 31 from the standpoint of the hydraulic power delivered 
and the brake efficiency. It shows that the hydraulic outpuc 
Stirling engine is capable of delivering 1.87 kW at 29.7 percent 
overall efficiency without the dynamic balance mechanism. 
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FIGURE 41. - Piston-displacer phase plot for the diaphragm piston 
plot displacer, hydraulic output Stirling engine 
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TABLE 11. - PERFORMANCE PREDICTIONS* OF DIAPHRAGM PISTON, 
FREE DISPLACER, STIRLING ENGINE WITH NULL 
CENTER BAND HYDRAULIC PUMP OUTPUT 



With the high phase 
displacer rod 
(Run 32d) 

With the new 
recommended 
displacer rod 
(Run 33) 

Displacer phase 
angle (deg) 

98 

77 

Indicated 
power (W) 

1210 

2010 

Indicated 

efficiency 

(percent) 

24.2 

32.0 

Power loss due 
to hydraulic 
system (W)** 

212 

212 

(due to dynamic 
balance 

mechanism only) 

(70) 

(70) 

Brake power (W) 



With balance 
mechanism 

998 

1798 

Without balance 
mechanism 

(1068) 

(1868) 

Brake efficiency 
(percent) 



With dynamic 
balance 

20 

28.6 

Without dynamic 
balance 

21.4 

29.7 

* Sunpower third order analysis, unconstrained 
transfer correlation and loss estimates. 
**See subsection 3.2.3. 

mode , FMA heat 
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Considering the fact that this engine design is not fully 
optimized, the level of performance as predicted here looks 
promising. 


3.2.3 Es timates of power losses due to the diaphragm assemblies 
and the hydraulic^ system . - The conversion oi RF-TOOO FPSE shown 
in Figure 1, into a hydraulic output Stirling engine as shown 
in Figure 18 introduced additional power losses for the engine. 
These losses are divided into three groups: 

« Pressure drop loss in various gaps, passages and 
valves . 

• Friction loss in the hydraulic side associated with 
the reciprocating components. 

• Leakage loss in the hydraulic system. 

These losses will decrease the hydraulic power output 
available for use and lower the engine brake efficiency. Each 
loss will be discussed in the following subsections. 


3. 2. 3.1 Pressure drop loss. - Table 12 lists the locations 
where the pressure drop loss occurs. The methods of calculating 
the pressure drop losses are briefly described below. Basically 
the analysis utilizes steady flow pressure drop correlations 
at each instant and location. 


For the gap between the diaphragm and its support, and 
for the inner cooler plug in the compression space, the in- 
stantaneous pressure drop, AP(t), between two points, Xj and 
X 2 # is obtai’ned from 



(4) 


where 

f 

P 

v 


d h (t) 


friction factor 
fluid density 

instantaneous, local fluid velocity 
instantaneous, local hydraulic diameter 
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TABLE 12. - ADDITIONAL PRESSURE DROP LOSS LOCATIONS 


Gas Side 

Inner cooler in the compression space 

The gap between the diaphragm and the gas side 
support (engine diaphragm and balance 
diaphragm) 

Hydraulic Side 

The gap between the diaphragm and the oil side 
support (engine diaphragm and balance 
diaphragm) 

Inlet valves of the pump 
Discharge valve of the pump 


For the valves, the instantaneous pressure drop is obtained from 
a standard orifice equation, such as 


AP(t) = £^4 
2 c d 

where 

C d = flow coefficient of a value 
p = fluid density. 

The loss of power, W, due to the pressure drops expressed by 
Equations (4) and (5) is obtained by 
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w 


1 

T 



AP ( t) q(t) dt 


( 6 ) 


where 

T = cycle period 

q(t) = instantaneous volume flow rate of fluid 

The integrands in Equations 4 and 6 contain the instantaneous 
local values for the velocity, the hydraulic diameter, and the 
flow rate. For a fixed geometry flow passage, such as the inner 
cooler, with known, d^, and known instantaneous mass flow rate, 
the velocity can be easily calculated from 


v = 


m 

PA 


(7) 


where 


m = instantaneous mass flow rate 
A = flow cross-sectional area 

By assuming that the valves of the hydraulic pump open and 
close instantly, the valves can also be viewed as fixed geometry 
flow passages for the duration of the valve flow. 

However, if the flow passage changes its shape constantly 
with time as is the case with the diaphragm gap, the expressions 
for v, d n , and q(t) become quite involved. For the diaphragm 
gap shown in Figure 43 the variables v, d^, and q are ex- 
pressed by 

v * c ll + C Sin ct) x [ x 2 2 + x ? x - 2x2 ] (8) 

2 h (x_ - x) 

d. = 2 — -£ (1 + £ sin cot) (9) 

n x « “ x * 


2 2 

(x 2 + x^x 2 - 2 x 1 ) cos cut ( 10 ) 


q = 


r h 

£ CO TT O 
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Figure 43. - Diaphragm support gap schematic. 
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where 


diaphragm design volume amplitude 
maxi mum theoretical volume ampl itude 

a) = rotational speed 

h Q = maximum diaphragm center deflection 
x = radial distance. 


Now the n? tantaneous pressure drop through the diaphragm 
gap, P (t) , end the resulting loss of power can be obtained 
from Equations 3 and 5 after substituting the expressions for 
V, dft, and q of Equations 8, 9, and 10. The remainder of the 
analysis is a straightforward mathematical procedure, and there- 
fore not given here. The results are instead given in Table 13, 
and will be briefly discussed. 

The high pressure drop loss of the inner cooler is a result 
of designing the inner cooler passage dimension to be smaller 
than the thermal diffusion layer thickness. The thermal diffusion 
layer thickness was calculated from the FMA cyclic heat transfer 
analysis mentioned earlier. 

The pressure drop losses in the diaphragm gaps of the gas 
side and the hydraulic side are negligible. The gaps were 
designed to have almost negligible pressure build up due to 
opening or closing of the gap. The pressure build up was cal- 
culated to be 690 Pa (0.1 psi) maximum for the two diaphragms. 

It was necessary to maintain a small pressure drop in order to 
prevent local deformations of the thin diaphragms. 


3. 2. 3. 2 Viscous friction and leakage loss. - The viscous friction 
loss means the hydrodynamic viscous friction occurring at the 
sliding surfaces of reciprocating components in the hydraulic 
system. The leakages occur past the same sliding surfaces. The 
locations for these losses are: 

• Transformer piston and rod. 

• Dynamic balance mechanism. 

• Null center band hydraulic pump and pump rod. 

Figure 44 is a schematic showing these locations which are 
indicated by thickened lines. The Roman numerals I to VI in- 
dicate the individual spaces that have different dynamic pressures. 
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TABLE 13. - PRESSURE DROP LOSS INTRODUCED BY THE 
CONVERSION OF RE- 1000 FPSE INTO 
HYDRAULIC OUTPUT ENGINE 


Locations 

Pressure drop loss 
(W) 

Gas side 

Inner cooler 

22 

Diaphragm gap 

Negligible (<1.0) 

Hydraulic side 

Diaphragm gap 

Negligible (<1.0) 

Slots in the support 

1.5 

Inlet valve 

10 

Discharge valve 

23 

Total 

. _ _ 

58 
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Figure 44. — Schematic of hydraulic system showing locations of viscous 
and leakage loss. 
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The calculation method of the losses at these locations is 
based on a steady laminar flow between stationary and moving 
flat, parallel plates. The shear force experienced by the 
plate moving at velocity V in x direction, at the distance 
b away from the parallel stationary plate, is given by 


Shear force = F = - ^ - P c 
s 2 dx M b 


( 11 ) 


where 


^ = pressure gradient in the flow direction 
y = viscosity 

The leakage rate, through the gap of width b, in the flow 
direction is given by 

Leakage rate =q=-^|£ + |v ,12) 


The power losses due to the viscous fricti n and the leakages 
are calculated from the following equations: 


Viscous loss 


W 


1 

Viscous T 



(13) 


Leakage loss = W 


Leakage 



AP dt 


(14) 


The integrands of Equations 13 and 14 include the variables 
F s , V, q, and AP, all of which are functions of time for each 
location shown in Figure 44. The reciprocating velocity V is 
given by the assumed Stirling engine dynamics: 


V = x n a) sin wt 


( 15 ) 


where x p = the displacement amplitude of the output piston. 
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At each location, the pressure aradient dP/dx is re- 
quired to calculate the two variables F s and q from Equations 
11 and 12. The dynamic pressures of spaces IV and V shown in 
Figure 44 were obtained by applying a dynamic force balance for 
each reciprocating components. The detail of the dynamic force 
balance will be omitted here. The resulting pressures are 
approximately given by the following equations. 


P_ ~ P 
I _ o 

P I1' P I1I — P H ! when Piping 

O : when not pumping 


x n 2 

P_„ c* P + - ° <*» s i n 

IV o 2 A_ 


P.. a; P_ - 


V ** *o 2 A 


m n X 1 

B ° sin u>t 


(16) 


P VI “ p o 


P VII “ p o + p c sin <ut " V 


where 

P o = mean engine gas pressure 

P u = high pressure hydraulic accumulator pressure 
m 0 * mass of the balance piston 

•* cross-sectional area of the balance piston 
P c = engine compression space pressure amplitude. 

Now, the pressures of the spaces I through VII given 
by Equation 16 can be used to calculate the shear force, F s , 
and the leakage rate, q, from Equations 11 and 12, respec- 
tively. Then, the shear force, F s , and the leakage rate, q, 
are substituted into Equation 13 to calculate the viscous loss 
and the leakage loss. The remaining detail will be again 
omitted and the results are given in T'.ble 14. 
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TABLE 14. - 


VISCOUS FLICTION LOSS AND LEAKAGE LOSS INTRODUCED 
BY THE CONVERSION OF RE-1000 FPSE INTO HYDRAULIC 
OUTPUT ENGINE 


Location 

Viscous friction loss 
(W) 

Leakage loss 
(W) 

Transformer piston 
and rod 

17 

5 

Dynamic balance 
piston 

58 

12 

Hydraulic pump 
piston 

17 

22 

Hydraulic pump 
rod 

3 

20 

Total 

95 

59 


Adding all three losses as shown in Table 13 and 14 the 
total additional loss caused by the conversion of the RE-1000 
FPSE into a hydraulic output engine is 212W. The optional 
dynamic balance mechanism accounts for approximately 70W of 
this total additional loss. 





















4. 


SUMMARY AND RECOMMENDATIONS 




4.1 Summary 

By modifying the Sunpower-buil t RE-1000 Free-Piston, Free- 
Displacer Stirling engine, the Foster-Mil ler/SUNPOWER team designed 
a hydraulic output version of the RE-1000 FPSE. The new engine has 
a diaphragm piston which is hydraulically coupled to a novel, null 
center band hydraulic pump as the output method. The pump pro- 
vides the inherent dynamic stability to the new engine. 

The engine is a test-bed engine rather than a prototype 
engine. Some of the notable features are listed below: 

• It is relatively simple in its design and inexpensive 
to build. 

• The modular design approach was used to facilitate the 
modification and replacement of the components. 

• A complete inertial balance of the moving parts is 
achievable in a horizontal or a vertical set up. 

• The start/stop operation is easily controllable. 

As a test-bed engine, the engine will be quite suitable for 
obtaining the following: 

• Cylinder heat transfer data in the narrow diaphragm 
space . 

• High cyclic rate diaphragm data. 

From the detailed analysis, the following observations were 

made : 

• For the simulation of the Free Displacer/Diaphraqm 
Piston Stirling engine with the null center band hydrau- 
lic pump output, the dynamically unconstrained third 
order Stirling engine analysis should be used. 

• The two critical parameters influencing the overall 
engine performance are the heat transfer in the dia- 
phragm space and the displacer phase angle. The heat 
transfer can be marginally controlled. However, the 
displacer angle can be changed fairly easily. 
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• In order to achieve a reasonable performance, the dia- 
phragm space should be actively cooled utilizing 
coolant passages, fins and flow patterns. 

• The proposed engine design is not optimized. To opti- 
mize properly, more knowledge in fluid flow and heat 
transfer phenomena near the diaphragm is required. 

Also, the optimization should include the entire engine: 
expansion space, regenerator, heater, cooler, diaphragm 
space, hydraulics, etc. 

• The diaphragm will last an indefinite period of time 
due to a very conservative design. The diaphragm edge 
supporting method using an ordinary rubber O-ring not 
only provides an effective sealing method but also re- 
duces the diaphragm stress for a given diaphragm volume 
displacement . 

Without the dynamic mass balancing mechanism, the engine will 
put out 1.87 kW of hydraulic power at 29.7 percent brake efficiency 
with a new, recommended displacer rod, and 1.07 kW at 21.4 percent 
with an existing displacer rod. The hydraulic output Stirling 
engine design concept can be easily scaled up to, say, 50 kW. For 
more power, any number of cylinders, can be ganged up. The hydrau- 
lic output is suitable for typical hydraulic accumulator powered 
systems at 100 to 200 bar range. 

4.2 Recommendations 

As a result of the design study, the following is recommended: 

• Build the hydraulic output Stirling engine with the new 
revised displacer rod design. This minor modification 
of the hot section will drastically increase the power 
and the efficiency of the engine. 

• Refine centering control schemes. 

• Test the hydraulic pump separately first. Modify and 
debug the pump before testing it in conjunction with 
the engine. 

• In the test plan, pay particular attention to the 
following : 

1. Diaphragm dynamics and positioning. 

2. Thermodynamics in the diaphragm region. 
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3. Displacer dynamics and the phase angle. 

4. Losses in hydraulic system. 

By varying the diaphragm midstroke position, it will be pos- 
sible to see the effect of the dead volume change on engine per- 
formance, especially since varying the rate and the temperature 
of the coolant flow in the diaphragm region will show the effect 
of heat transfer on the engine performance. 


125 


5. REFERENCES 


Martini, W.R., Stirling Engine Newsletter, Martini Engi- 
neering, Richland, Washington, November, 1981. 

Kalnins, A., and Lestingi, J.F., "On Nonlinear Analysis of 
Elastic Shells of Revolution," ASME Trans., Vol . 34, Series 
E, No. 1, pp. 59-64, March, 1967. 

Wahl, A.M., Mechanical Springs, McGraw-Hill Book Co., 

Second Edition, pp. 170-172. 

Allegheny Ludlum Data Sheet 115-92358-350, Allegheny Ludlum 
Steel Corporation, Research Laboratory, Brackenridge , 
Pennsylvania . 

Taylor, C.F., The Internal Combustion Engine in Theory and 
Practice, Second Edition, MIT Press, Cambridge, Massachusetts 
pp. 266-311. 

Lee, K., Smith, J.L. Jr., and Faulkner, H.B., "Performance 
Loss Due to Transient Heat Transfer in the Cylinders of 
Stirling Engines," Proceedings of the Fifteenth IECEC 
Conference, Seattle, Washington, August, 1980. 

Foster-Miller Associates, Inc., "Design and Development of 
Stirling Engines for Stationary Power Generation Applications 
in the 500 to 3000 hp Range, " Phase I Final Report to 
Argonne National Laboratory, Argonne, IL, No. DE-AC02- 
79ET15208, October, 1980. 


EBECED1NG PAGE BLANK NOT FILMED 


127 


